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Pelos inícios do Novo Ano, 

surgiu o primeiro número da Técnica de 1996, 
sinal de vida, 

de vida sincronizada com o tempo. 

Em 1996 completa o Instituto Superior Técnico 85 anos. | 
Sob a orientação do seu Presidente, 

uma comissão diversificada 

trabalha 

no sentido de uma justa e condigna 
comemoração dessas efemérides. 

Gentilmente, 

a Técnica foi convidada 

a comparticipar. 

Por ocasião dos 75 Anos, 

num número especialmente significativo, 
debruçou-se a Técnica 

sobre os antecedentes, origens, 

no que de importante e nesse período decorreu 
no Instituto. | 
Pela ocasião do 80 anos do IST, prospectou-se 
o futuro 

das Ciências e Tecnologias 

de que esta Escola se ocupa e ocupará. 
E. mais uma vez, 

a Técnica foi chamada a comparticipar. 
Hoje. 

por ocasião dos 85 Anos, 
sente-se que o determinante 

e característico, 

é O que. organizada 

e sistematicamente, se investiga 

e quem nisso se investiga 

e quem nisso se envolve. 

Uma notícia sucinta e rigorosa 

é o que se pode e quanto bastará, 
para construir 

não exactamente uma fotografia, 
mas uma mancha, 

suficientemente indiciadora, 

da Investigação e Desenvolvimento 
no Instituto Superior Técnico 

por ocasião dos seus 85 Anos. 


| 
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SUMARIO 


Pulsares são estrelas relativis- 
tas, cujos mecanismos de radi- 
ação se suportam nos modelos 
das respectivas magnetosferas. 
Embora a estrela seja com- 
pacta e a magnetosfera exija 
uma descrição tendo em conta 
todos os aspectos geometrodi- 
nâmicos, é possível isolar os 
efeitos puramente electromag- 
néticos, no contexto do espaço 
tetradimensional de Minkowski. 
À aproximação de particulas 
singulares do plasma magne- 
tosférico é seguida, e efeitos 
radiativos têm de ser conside- 
rados na equação do movi- 
mento. Então uma descrição 
quasi-clássica é necessária, 
com trajectórias deterministicas 
substituídas por trajectórias 
estocasticas, pois correcções 
quânticas devidas aos efeitos 
radiativos têm de ser conside- 
radas. Porém ao considerar os 
efeitos de Relatividade Geral o 
movimento do plasma será 
diferente. Neste artigo serão 
discutidos os processos físicos 
associados. 


ABSTRACT 


Pulsars are relativistic stars, 
whose radiation mechanisms 
are related to the models of 
their magnetospheres. Although 
the star is compact and its 
magnetosphere requires a des- 
cription taking into account all 
geometrodynimical aspects, it is 
possible to isolate the pure 
electromagnetic effects, in the 
context of Minkowski 4-space. 
The single particle approach is 
followed, and radiative effects 
must be taken into account in 
the equation of motion. Then a 
quasi-classical description is 
required with deterministic tra- 
jectories substituted by sto- 
chastic ones, as radiative quan- 
tum corrections have to be 
considered. Nevertheless when 


GEOMETRODINÂMICA 
DE PULSARES: 
AS MAGNETOSFERAS 


A.A. DA COSTA 
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INTRODUÇÃO 


A descoberta em 1967 por Jocelyn Bell (Hewish et al, 1968) do objecto 
PSR1919+21 criou um novo tipo de objectos cósmicos, os Pulsares, que se 
manifestam pela emissão de radiação electromagnética não térmica, modu- 
lada em impulsos periódicos. Formam dois grupos: o Grupo |, com períodos 
0.033s < P < 3.05 (o PSR1919+21 tem P = 1.335); o Grupo Il apresenta 
períodos ainda mais curtos, 1.55ms < P < 10ms (Camilo, 1994). Existem 
ainda em sistemas binários, tanto do tipo | como do tipo Il e radiam funda- 
mentalmente na banda de radio frequência. Porém os pulsares do Caran- 
guejo (CP), P=0.033s, ou de Vela, P=0.088s entre outros foram detectados 
noutras bandas electromagnéticas, e.g. óptica e raios-y. 


Estes impulsos individuais são morfologicamente muito diferentes uns dos 
outros duma forma aleatoria, mas com uma grande regularidade estatística, 
pois a média de cerca de 500 impulsos em sucessão no tempo, apresenta 
sempre o mesmo impulso integrado qualquer que seja a amostra (Fig. 1). 
Os impulsos individuais apresentam ainda elevado grau de polarização que 
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Fig. 1 Registo de 100 impulsos de PSR1133+66 a 600MHz — Impulsos consecutivos são 
mostrados verticalmente para mostrar variações de forma e tempos de chegada. Uma 
média de 500 impulsos, o impulso Integrado, é mostrado no topo, e é fixo (segundo Cor- 
des, 1979, por cortesia do autor). 


Eai) 


diminui com a frequência, embora seja uma polarização modulada, isto é, 
varia ao longo do impulso, e de impulso para impulso, verificando-se que a 
polarização do impulso integrado é inferior à dos impulsos individuais. Nal- 
guns pulsares com uma resolução temporal da ordem dos microssegundos, 
os impulsos desdobram-se em microimpulsos (Fig. 2). 


FULL 
BANDWIDTH 


to mem 
IO msec TIME 


Pd W 318 MHz 
dr 
[zo fu ga” 


is a PORN 


[so fu ta id TOPA, HH MHz 
a a Ms 
ci 5 msec ta 


DT] 


Fig. 2 Microestrutura e impulsos integrados - No topo mostra-se um impulso do 
PSR0950+08 observado a 111.5 MHz. O tempo de resolução é de 56 ms. Porém o Impulso 
integrado formado a partir de 1500 impulsos individuais, é semelhante a estes, quando 
estes são observados com tempos de resolução normais da ordem dos milissegundos. 
Neste caso é interessante observar o “ruido” sobreposto ao Impulso Integrado consequên- 
cia da microestrutura (segundo Rickett et al, 1975; por cortesia dos autores) 


A descoberta do CP, mostrou que os pulsares eram estrelas de neutrões 
em rotação rápida, de acordo com a Hipótese de Gold (1968). Só estes cor- 
pos compactos com uma massa inferior a 3M; (em unidades geométricas 
1Mo- 3 km) e um raio da ordem dos 10 km, permitiam explicar não só os 
curtos períodos observados, mas também a sua enorme variedade e preci- 
são. São os relógios mais precisos. 


Admite-se que as estrelas de neutrões do Grupo | são o que resta das 
supernovas pela implosão do núcleo das estrelas que explodem. A implosão 
conserva o momento angular o que implica um período inicial P <P, 
devido ao processo posterior de desaceleração da estrela; e conserva o 
fluxo magnético resultando uma indução magnética final à superfície B - 


10ºT, valor assumido e associado a um momento dipolar magnético d. Ha- 


veria assim um campo electromagnético (CEM) girante e dipolar corrodando 
com a estrela com velocidade angular S2, e, por isso mesmo, de baixa fre- 
quência, o modelo de Deutsch generalizado (Deutsch, 1955), desde que o 
espaço à volta da estrela fosse vácuo electromagnético, e o espaço-tempo 
associado, o espaço quadridimensional de Minkowski (espaço-M). A potên- 
cia radiada, com B, - 10ºT, é igual à potência cinética de desaceleração 


1202, onde 1 é o momento de inércia e Q tem um valor proveniente da 
observação. Porém o mecanismo de manutenção deste CEM é desconhe- 
cido. Admite-se que a indução magnética enfraquece ao mesmo tempo que 
P cresce, e para P > 35 tem um valor abaixo dum valor crítico para o qual 
os pulsares deixam de radiar. As respectivas estrelas de neutrões passa- 


General Relativistic effects are 
considered, the motion of 
plasma will be different. In this 
paper the associated physical 
processes are discussed. 


riam por um largo período mortas e ressuscitariam nos pulsares do Grupo | 
ao receberem uma injecção de matéria e respectivo momento angular, e 
voltavam a radiar com os períodos muito curtos observados e 5B, = 10“T 
(para mais detalhes ver Dewey, 1989). 


O valor de 5 define o exterior da estrela como uma região extensa, a mag- 
netosfera. Esta região não é vácuo electromagnético, mas contem um 
plasma de electrões e positrões cuja fonte é a estrela, embebido num CEM 
global girante de baixa frequência £2, sobreposição do CEM dipolar com um 
CEM não dipolar de que o CEM dipolar é raiz, gerado pelo movimento do 
plasma magnetosférico no CEM global. A determinação deste é assim um 
problema funcional complexo, pois depende da relação mútua entre a estru- 
tura do CEM que acelera o plasma e a distribuição deste, e da distribuição 
do CEM a superficie. 

Este plasma existe devido à condição fronteira electromagnética à superfície 
da estrela. Simplificando, no interior manifesta-se a condição de condutivi- 
dade perfeita, que em espaço-M é E+(Qxr)xB=0; mas no exterior as par- 
tículas podem mover-se ao longo das linhas de força do campo magnético e 
assim E+(Oxr)xB=-Vy (Goldreich & Julian, 1969; Mestel, 1971). Existirá 
uma densidade de carga à superfície da estrela, o campo eléctrico terá aí 
descontinuidades devido à rotação, e como E.B=() à superfície, cargas são 
retirados da estrela, processo físico que não está claramente estabelecido 
(Ruderman & Sutherland, 1975), e lançados na magnetosfera. Forma-se 
assim um plasma de electrões positrões separaveis naturalmente sublumi- 
nosos, com uma densidade de corrente de condução e carga de espaço 
associadas, rodando com a velocidade angular do pulsar, e que quando 
situadas para além do cilindro de luz do pulsar (CLP), cilindro virtual que 
corrodando com a estrela terá a velocidade da luz, e raio R,=c/£2, apre- 
sentam velocidades lineares azimutais superluminosas, sendo isto possivel 
porque não possuem inércia. Estas distribuições geram a componente não- 
dipolar do CEM, devido não só ao valor de E.B/8+0 em toda a parte, que 
uma tal componente procura diminuir (Princípio de Le Chatelier, Mestel et al, 
1976) alterando também B, mas também devido ao comportamento diná- 
mico diverso das diferentes partículas do plasma, relacionado com a geome- 
tria do CEM global, e valor da relação E/cB. Estas afirmações são verdadei- 


ras mesmo para o caso em que Sid, e assim mesmo no caso axissimétrico 
a estrela de neutrões é activa. No caso geral 0 d+0, o CEM não dipolar 
não tem simetria azimutal e a longa distância do pulsar constitui a radiação 
de radio frequência. 

Desta forma pulsares são estrelas relativistas que contém fenomenologia 
que se desenvolve em condições extremas, pois são corpos hipercompac- 
tos, com R./2M<3 em unidades geométricas (2M=1M,;), com CEM hiperin- 
tensos. Assim o interior da estrela e a magnetosfera constituem dois espa- 
ços geometrodinâmicos (espaços-G), cujas métricas têm o mesmo valor à 
superficie, em particular no caso de pulsares do Grupo Il com um R,-4R.. 
A estrela apresenta fenomenologia diversa no seu interior, à sua superfície, 
e no seu exterior, não sendo suficientes descrições clássicas. As descrições 
quânticas em espaços-G são inevitáveis. 

O interior apresenta densidades internas muito elevadas no seu centro, 
superiores à densidade dos núcleos atómicos (supra-nucleares), limitada 
superiormente pela impossibilidade da velocidade do som não poder ser 
superior à velocidade da luz; mais próximo da superfície as densidades 
decrescem, tornando possivel matéria formada por neutrões e possuindo um 
electrão/protão por cada 8 (oito) neutrões. O interior da estrela apresenta 
assim superfluidez e supercondutividade. A superfície será uma camada de 
transição entre o interior e o exterior da estrela, definindo uma condição 


fronteira difusa, com uma crosta formada por Fe* e encontra-se no estado 
sólido apesar da sua temperatura ser da ordem dos 10ºkK. Como B - 
015,,0onde B -4xI0'T é limite superior da aproximação classica então 
a superfície é constituida por uma rede atómica polimerizada de atomos que 
em campos magnéticos intensos, apresentam uma nuvem electrónica cilin- 
drica, e que devido a esta estrutura, encaixam uns nos outros, permitindo a 
existência de camadas de electrões livres que podem ser facilmente extraí- 
dos (Irvine, 1978). O movimento próximo da superfície da estrela tem um 
limite superior de energia, devido ao efeito superluminoso de Ginzburg 
(Ginzburg, 1989) generalizado a espaços-G. Ele foi estabelecido para espa- 
cos-M relacionado com polarização do vácuo, que começa a ser importante 
quando y5B>1I0'T;, transformando o vácuo num meio anisotrópico, e que 
impede valores elevados de y (o factor de Lorentz), necessários para expli- 
car a radiação-y. 

O plasma magnetosférico, constituído por electrões e positrões, devido ao 
CEM hiperintenso apresentará regiões ultrarelativistas, cuja distribuição tem 
em conta as condições fronteiras electromagnéticas e geometrodinâmicas. 
Neste trabalho far-se-á o estudo das magnetosferas na vizinhança do CLP, 
para pulsares do Grupo |, em espaço-M, pois o CLP tem um raio directa- 
mente proporcional a P, com 1500 km para o CP. Assim usar-se-à electrodi- 
nâmica clássica com correcções quânticas, admitindo condições fronteiras à 
superfície compatíveis. Concluir-se-á que a radiação no domínios das radio 
frequências e dos raios-y poderá existir na vizinhança ou no exterior do 
CLP. Será ainda feita uma generalização a espaços-G, identificando tam- 
bém os problemas associados a uma tal descrição. 


2. A DESCRIÇÃO MINKOWSKIANA DAS MAGNETOSFERAS 


2.1. Movimento de Partículas numa geometria bi-dimensional 


Num pulsar d.2 70, combinação não-linear de: d/O, o modelo axissimé- 
trico; e di£?, o modelo perpendicular. Neste artigo situar-nos-emos no 
modelo perpendicular, o modelo mais simples que permite explicar a radia- 
ção observada, por não possuir simetria azimutal. O modelo axissimétrico 
pode ser consultado em Fitzpatrick & Mestel (1988,a,b). 


O CEM pode ser assumido como dipolar, desde que consideremos um 
plasma de baixa densidade. Então o plasma é um conjunto de partículas 
isoladas, consequência da energia electromagnética do CEM ser muito 
superior à energia de interacção das partículas, o que reforça a utilização da 
Aproximação de Partículas Singulares (APS). Porém fora do plano equa- 
torial o CEM tridimensional dipolar tem E.Bz0. Simplificando (da Costa & 
Kahn, 1982; da Costa, 1983), usa-se uma configuração bidimensional. Esta 
estrutura cilíndrica parecer-se-á com a configuração tridimensional no plano 
equatorial, e os seus resultados farão sentido se nos mantivermos numa 
regiao pequena, com uma dimensão linear da ordem de 0.1&R,. 


Numa geometria 2-D o CEM “dipolar' é derivado do potencial vector 
r(B) 


A=Aydi (B)+Y;(B) cos PD + ar nd e a (1) 


Nesta expressão J, e Y, são funções de Bessel, a |2, B=p/R e PY é 


um parâmetro de fase. Com este potencial vector o CEM tem EJA, HIA e 


as linhas de campo magnético (LCM) são linhas de A constante A =A,C., 


ss + SEO E CT TT OO o e 


que define P=qp-wr='P(5). O valor de 'P, é escolhido de maneira que 


C.= J(1)+Y;(1)=0.8966 para a LCM tangente ao CLP, e W'=0 para 


P=0 e B=1, que corresponde à LCM C.=0.6858 (fig. 3). 


Fig. 3 As linhas de força do campo magnético do modelo bidimensional —- Nesta ima- 
gem estão representadas as regiões E<cB, E>cB, as regiões de transição E=cB entre uma 
e outra e o segmento do CLP onde as linhas de força se orientam na direcção do movi- 
mento. Distingue-se ainda a região S, , ea região S.' (ver texto). O rectângulo inferior 
mostra com mais detalhe a região do CLP onde as linhas de força se orientam na direcção 
do movimento. A cheio estão as LCM Cc, = 0.66, 0.6858, 0.70, 0.75, 0.BO. C. =0.6858 é 
perpendicular ao CLP, C, =0.8966 é tangente ao CLP (segundo da Costa & Kahn, 1982; 
e da Costa, 1983). 


As partículas iniciam o seu movimento junto à superfície da estrela com : - 
0 e y,=1, próximo de uma LCM. A partícula move-se ao longo das LCM. 


Porém quando as partículas atingem a secção do CLP aonde W'>0, ou 

quando atingem as duas superfícies de transição das regiões 

E<cB>E>cB, superfícies S:> e S-, a energia e>00 e a partícula 

tem de abandonar a LCM a muito curta distância das superfícies, da ordem 

dos 10º K para o CP, Este valor tão baixo é consequência do parâmetro 
EA, 


(p= 
RC (2) 


que mede a intensidade do campo electromagnético (Rees, 1971), razão 
entre a energia adquirida por uma partícula sujeita a uma diferença de 
potencial numa distância da ordem de R,,e a energia de repouso da 


mesma. que é inversamente proporcional a P? e para o CP e electrões vale 
00. a 

n=10" junto ao CLP onde B - 100T. Um cálculo tem então de começar, 

que tem de ter em conta a reacção de radiação sobre o movimento das par- 

tículas. 


2.2. A Aproximação Quase-Clássica 


Os estudos sobre os efeitos da radiação das partículas sobre a sua própria 
trajectória remontam a Schott (1912), mas a sua formulação covariante no 
domínio clássico foi feita por Dirac (1938). A descrição clássica é válida 
desde que ho, <<ymc, com w,..=y'c/p, onde r é o raio da trajectória, 
o valor de w para o qual o espectro da potência radiada P,., expresso por 
Poslp, n, w, t), é máximo. A equação de Lorentz-Dirac do movimento 
(EL-DM) duma partícula clássica é 


W = Fu, Exu + A (uu Je 
m ÇA (3) 


com c — 6.27 x 10“'s para um electrão (positrão), o tempo que um sinal leva 
a percorrer uma distância igual a 2/3 do raio do electrão/positrão, e onde os 
três termos na EL-DM são respectivamente os termos de Lorentz, Schott e 
Dirac. O termo correspondente à coordenada tempo mostra a possibilidade 
de uma partícula radiar tanta P,., como aquela que adquire através da força 
de Lorentz, P,... Nesta situação a partícula tem y - 10” e radia na banda dos 
raios-y duros, da ordem dos 10GeV. A descrição do campo radiado têm de 
ser feita por fotões, a energia é radiada de forma discreta e aleatória, e dai 
que os seus efeitos são também discretos e aleatórios. Porém a EL-DM foi 
estabelecida para campos radiados com descrição ondulatória, sendo a 
energia radiada de forma continua. Logo não pode ser usada. 

Nestas circunstâncias tem de ser usado o regime q-classico (quase-clássico), 
uma correcção quântica ao regime clássico (da Costa & Kahn, 1991). Nesta 
descrição tal como na correspondente descrição clássica sem efeitos radiati- 
vos, as partículas movem-se devido à acção da força de Lorentz. Porém há 
mudanças bruscas de momento e energia devido à radiação de fotões de 
raios-y com frequências aleatórias, a intervalos de tempo aleatórios. As tra- 
jectórias passam de deterministicas a estocásticas. Trata-se de um processo 
de Poisson composto, com A(p. ?) a probabilidade de radiar um fotão por uni- 
dade de tempo, e que é igual ao número médio de fotões radiados por uni- 
dade de tempo. Estes intervalos de tempo t obedecem à distribuição 


FlT)=A(pai+rT) exp/ [di A(p«)] (4) 


Estes fotões individuais têm uma distribuição de densidade de probabilidade 
não normalizada P (p.n.orj=(h0) P (p.n.o.), com momento linear p, e 
radiando kk=h(w/c)n, no interior dum cone cujo eixo é paralelo a p, e a 
abertura y''. p é assim o eixo polar dum sistema de coordenadas esféricas. 
Desta forma k define três variáveis aleatórias: w, 0 e q. 


A radiação de fotões individuais estabelece uma fronteira difusa entre os 
domínios clássico e q-clássico onde how-mc: esta expressão apresenta 
uma extraordinária semelhança com a fronteira entre regimes não-relativis- 
tas e relativistas em plasmas térmicos onde &T -mc” . Ao valor acima de w 
chamaremos q, ,,(- 10 Hz) e daí que fotões com v<w,,, não afectam o 
movimento das partículas. As partículas radiam preferencialmente em Wax 
e fotões com w>>ko,.(k>1), têm uma baixa probabilidade de ocorrerem. 
Portanto só quando O, <Q, há uma probabilidade da trajectória da par- 
ticula ser afectada pela radiação. 

Este processo de Poisson composto é de facto também um Processo de 
Markov discreto e é um caso particular dum processo de pontos com mar- 


qu a 


cas (Snyder & Miller, 1991), desde que tw zmc . Definamos o espaço de 
marcas | como o espaço de ocorrências de (mw, 0, q). Seja To intervalo de 
tempo [f, =))e 7xu o espaço onde todas as marcas ocorrem. Então M (T, 
B) é o número de pontos existente no espaço TxB, onde Tcf, e nesta 
nova descrição q-clássica a equação do movimento é 


dp(t)= e(Elr(t).t)+ v(r)x B[r(1).]jdr — | tkfr(o) pl: ). 8, U]M(dr,dU) (5) 


A diferencial dp(1)=p(t+dr)-p(!) é o incremento em p(t) que ocorre 
durante (t, tr+dtf]. Colocando 


(0, nenhum fotão radiado no tempo dt, ou, se radiado 0 < 0, 


| m(ar,dU)= un 
+ lt. o fotão radiado no tempo dr tiver 0 > 


e (6) 


Se nenhum fotão é radiado, ou se radiado, v<o,,, então 
dp(r)=e(E+vxB)dr; mas se um fotão é radiado com frequência 0>0,.,, 
então um ponto com uma marca significante u ocorre, e 
dp(r)=e(E+vxB)dr-kk[r(1).p(r).r.u|. Como só há descontinuidades nos 
pontos de radiação eficazes, a integração é feita de ponto a ponto de radia- 
ção eficaz, só com a Força de Lorentz, e no extremo do intervalo deduz-se 
o momento linear do fotão radiado. 


Para a magnetosfera do CP, a fronteira entre domínios clássicos e q-clássi- 
cos é atingida quando y - 10º. A simulação na região q-clássica suporta-se 
num Método de Monte-Carlo para gerar as várias variáveis aleatórias em 
presença. No caso dos pulsares interessa cálcular as pseudo-trajectórias no 
referencial que roda com a estrela, mas com as grandezas de momento e 
energia no referencial do observador. Os resultados (Fig. 4) mostram que 
várias particulas, com a mesma pseudo trajectória clássica, seguirão pseudo 
trajectórias aleatórias todas diferentes, mas muito próximas umas das 
outras. Isto sucede para todas as particulas, e portanto é expectável um 
feixe de trajectórias atingindo y - 10”. Para esse valor P,.,/P,, -1. À tempe- 
ratura de um plasma equivalente com a distribuição de energias apresen- 
tada por este feixe de pseudo-trajectórias seria cerca de T,, -I0ºK, o que 
mostra que o processo radiativo é largamente dominante, o que reforça a 
APS, pois as interacções entre partículas não são importantes. Quando as 
trajectórias estão bastante afastadas da região onde E>cB(C. 20.75), há 
uma influência do movimento giromagnético nas trajectórias, não dominante, 
visto que a rotação do CEM de baixa frequência domina o processo (Kirk, 
1980). Neste caso particular há fenómenos de ressonância observados nal- 
gumas trajectórias, que amplifica a perturbação giromagnética. Todos estes 
fenómenos se passam para distâncias muito curtas, da ordem dos 
10“ =10* -10R,, a que correspondem distâncias da ordem dos 10 - 1 - 0.1m 


conforme as trajectórias e o valor de C., pois quanto maior C.>0.75, 
menor a escala do fenómeno. 
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Fig. 4 Comparação entre pseudo-trajectórias para (.=0.75 para o caso clássico e 
quase-clássico — No topo (figura A) mostra-se a representação do espaço aonde são fei- 
tas os cálculos das pseudo trajectórias para uma CLP que exista só em regiões onde 
E<cB. No meio (figura B) mostra-se o cálculo duma pseudo-trajectória usando EL-DM. Na 
ultima figura faz-se o mesmo cálculo para 6 pseudo-trajectórias usando considerações 
semiclássicas. A diferença é notória. Repare-se que enquanto para EL-DM todas as 6 
pseudo-trajectórias coincidem o mesmo não se passa para O caso q-clássico, embora 
exista colimação (segundo da Costa, 1995). 


2.3. O Comportamento Estatístico 


A descrição do comportamento do plasma enquanto entidade colectiva, 
suporta-se na determinação da função de distribuição cinética /(r,.p.!), 
onde r, contêm as coordenadas no referencial de rotação, e que deve ser 
considerada como descrevendo um comportamento médio do plasma radia- 
tivo da magnetosfera. Como um tal plasma tem uma T, -I0ºK, a APS 


será a base do cálculo de /lr,.p.!) (Kahn & da Costa, 1984), se tal descri- 
ção for possivel. 

O cálculo de exige frequências aleatórias e discretização no espaço. Como 
as trajectórias são estocásticas, esta discretização terá de ser inferior à 
dimensão minima do fenómeno estocástico na direcção do movimento, e 
muito menor que o espraiamento das trajectórias no plano normal à trajectó- 
ria, e nessas células deverá existir um número de partículas suficientemente 
significativo para que o plasma ainda possa ser identificado como tal. A den- 
sidade necessária à radiação de raios-y é da ordem de 10" — 10" electrões 
m (da Costa & Kahn, 1983). Por exemplo, para C.=0.75 tal implica a pre- 


sença nas células de 10'-I0" partículas. Neste caso existe um processo 
estocástico médio regular, e assim é possivel estabelecer variedades regu- 


lares médias, formadas por trajectórias. Porém quando C.>0.75, para cada 
ordem de grandeza que a dimensão linear da célula decresça, o número de 
partículas por célula diminui de três ordens de grandeza. O carácter colec- 
tivo do plasma enfraquece. 


Atendendo ao carácter descontínuo da radiação de fotões, haverá uma com- 
ponente aleatória que varia no tempo e que deverá ser amplificada com o 
movimento do plasma. Este efeito é tanto mais forte quanto mais fraca for a 
descrição colectiva do plasma. A APS implica trajectórias estocásticas e 
deste modo o cálculo de flr,.p.t) terá um carácter recursivo contínuo, com 
bifurcações múltiplas em toda a parte recursiva do cálculo. É assim possível 
encarar o plasma magnetosférico como um meio aleatório com característi- 
cas fractais. 


No modelo bi-dimensional os electrões e positrões têm comportamento 
semelhante, mas quando :=>c os electrões e positrões afastam-se do 
plano equatorial em direcções opostas. Forma-se uma distribuição de densi- 
dade de corrente que pode rodar com velocidade superluminosa se estiver 
além do CLP. Veremos posteriormente quais as consequências. 


2.4. O Modelo Tridimensional 


A teoria exposta é independente dos modelos, e aplicada depois ao caso 
bidimensional. Portanto no modelo 3-D com £2.d*0 todos os fenómenos 
assinalados deverão estar presentes, mas as distribuições de carga de 
espaço e de densidade de corrente de condução estarão presentes duma 
forma ainda mais intensa, devido à assimetria de comportamento dos elec- 
trões e positrões devido ao ângulo entre Q e d, e à presença de E.Bz0, 
Num tal modelo, ainda por desenvolver, as condições fronteiras à superfície 
terão um papel altamente selectivo. 


Sabe-se entretanto que as distribuições de densidades de corrente e de 
carga de espaço, sem inércia, ao estenderem-se para além do CLP, corro- 
dam com a estrela movendo-se com velocidades superluminosa. Em mag- 
netosferas de pulsares o movimento do ponto-fonte das distribuições de car- 
gas e correntes é circular e temos um número ímpar de ondas associadas a 
cada ponto-campo, consequência da expressão 


i (7) 


onde R é a distancia do ponto-campo ao ponto-fonte, poder ter mais de uma 
solução para cada t (fig. 5). Se o raio da circunferência da trajectória do 
ponto-fonte for B<1 ou o ponto-campo for interior ao CLP, então existe sem- 
pre uma e uma só onda no ponto-campo. Mas se B>1 e o ponto-campo esti- 
ver fora do CLP, então (7) dá dois números impares de ondas consecutivos, 
sendo o par função do raio da circunferência do ponto-fonte, estando o 
ponto-campo a receber um ou outro dos valores conforme o instante de 
recepção. A partir do CLP, com b crescente, o primeiro par é (1,3), depois 
muda para (2n-1, 2n+1) até = sendo n é a ordem do par tal que 
B E Bacia ), sendo os extremos do intervalo as raízes de ordem n e n+1 


da equação 
tg vB —| = VB —1. (8) 


O CEM gerado pela distribuição de carga de espaço e densidade de cor- 
rente é assim semelhante, mas diferente do efeito Cerenkov. A superfície 
singular associada a cada ponto fonte é a justaposição duma superfície de 
Cerenkov associado aos mínimos das curvas t=t(s), com uma superfície cor- 


Fig. 5 Quatro diferentes representações de t = s+R(s)/c para quatro diferentes valores 
do raio da trajectória - O caso para a existência de mais de uma onda pode ser testemu- 
nhado para b= 2.5 e b,=6. b,=1 é um caso degenerado (segundo da Costa & Kahn, 1985) 


a respondente aos máximos, a superfície de Schott: é a superfície de Schott- 
Cerenkov (fig. 6; da Costa & Kahn, 1985). Também neste caso um contínuo 
de pontos-fonte fará desaparecer as singularidades assinaladas. 
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Fig. 6 Uma volta da superfície de Schott-Cerenkov — A figura mostra a intersecção da 

superficie com planos a altura constante. A superfície de Schott está escondida pela 

superficie de Cerenkov e está representada por linhas a traço interrompido. Como se pode 

ver a superficie de Schott desenvolve-se a partir duma cúspide no plano de simetria. O cir- | 
culo a traço interrompido representa a trajectória do ponto-fonte cónico raio b=3.0 
(segundo da Costa & Kahn, 1985). 


O efeito de Schott-Cerenkov, um fenómeno superluminoso, para além de 
explicar as alterações não-dipolares ao CEM dipolar de baixa frequência do 
pulsar, deverá ainda explicar a radiação de radio-frequência (da Costa, 
1976; Ardavan, 1994). As distribuições de densidade de carga e corrente 
associadas variam em distâncias d<<R,, e por isso a radiação de radio fre- 
quência constitui uma generalização da teoria da radiação sincrotrónica. 
A única diferença está no facto de que a direcção da radiação sincrotrónica 
máxima estar no plano osculador da trajectória, e no caso de radiação 
superluminosa, estar para fora do mesmo plano. 


Como estas distribuições são também estocásticas, geram CEM de baixa 
frequência estocásticos, reforçando assim a estocacidade do processo. 
Por outro lado a radiação de radio frequência vem também afectada. 
Estes factos podem explicar as diferenças entre os impulsos individuais, e 
os efeitos de microestrutura (ruído) observado nalguns pulsares (Rickett et 
al 1975), e convergem com os trabalhos de Zhuravlev & Popov (1992), 
que identificaram componentes estocásticas no processo de radiação de 


radiofrequência. 


3. A DESCRIÇÃO GEOMETRODINÂMICA DA ESTRELA 


3.1. A Métrica da Estrela 


Numa descrição geometrodinâmica, a determinação das métrica é o pri- 
meiro problema a ser resolvido. A métrica fora da estrela depende da densi- 
dade da matéria no seu interior, do seu momento angular e da energia do 
CEM de baixa frequência associado à estrela. A métrica interna influencia a 
equação de estado, e portanto, a densidade da matéria dentro da estrela. 
Ambas as métricas têm o mesmo valor à superfície da estrela, e desta 
forma influenciam-se mutuamente. Estamos perante um problema funcional 
do qual não existem soluções exactas. 

Considerando só mecânica newtoniana, a energia cinética de rotação do 
pulsar do Caranguejo é TT. =0.5/0' - 10], e a energia electromagnética 
da estrela é W.,-10"]. Como W,.=GM/R-10“J, então a estrela seria 
esfericamente simétrica se 7. =0, e axialmente simétrica segundo OQ sse 
Old. Repare-se que 7, e W,, são correcções a W... Esta é para 1Ms 
cerca de 15% da massa bariónica da estrela, cuja massa real é assim 
menor do que a sua massa bariônica, sendo o que falta toda a energia per- 
dida no processo de formação do pulsar, devido à radiações luminosa, de 
neutrinos ou gravitacional (Kippenhahn & Weigert, 1990). A métrica fora da 
estrela não pode em geral ser a métrica de Kerr (Andrade, 1994), devido à 
existência da métrica interna, situação diversa dum buraco negro. Porém 
tanto 7, como W, são muito menores que a correcção devida ao campo 
gravítico. A métrica da estrela é por conseguinte uma métrica de Schwars- 
child perturbada (MSP), ou seja é uma série de potências de duas variá- 
veis, uma e=(9/0,. F | onde 92, é o valor de Q onde a força centrifuga 
iguala a força gravitacional no equador da estrela (Hartle, 1967); outra uma 
variável associada a uma energia electromagnética limite igual à energia 
gravitacional da estrela; o termo de ordem zero é a métrica de Schwars- 
child. Note-se que e<1 sempre; e para o cálculo do termo de ordem zero a 
massa a considerar é tão só a soma das componentes bariónica e gravita- 
cional. Se T,, é uma perturbação de segunda ordem, então W., é de 
quarta ordem. 


ti 


13 


14 


3.2. Problemas Associados com a Descrição Geometrodinâmica 


Numa descrição em espaços-G, o comportamento electrodinâmico do 
plasma magnetosférico depende do CEM de baixa frequência, da métrica 
externa e das condições fronteiras à superfície; a influência do interior da 
estrela tem de ser considerada. Os efeitos de rotação diferem consideravel- 
mente da teoria Newtoniana, e são (Miller, 1992): 


1. O arrastamento dos referenciais de inércia; 
2. Mudança nos efeitos centrifugos; 
3. Mudança da relação entre o momento angular e a velocidade angular. 


O tema 2, menos conhecido, tem sido estudado no contexto de uma métrica 
de Schwarschild com uma partícula rodando numa órbita circular à sua volta 
(Abramowicz, 1990; Abramowicz et al, 1993). Se uma densidade constante 
de matéria é assumida dentro do objecto, existe uma região aonde a força 
centrifuga se dirige na direcção do interior da estrela, para R/2M<l.5 (Mil- 
ler, 1992; Abramowicz, 1990). O mesmo se espera que aconteça para pul- 
sares, pois a métrica é a MSP e embora aqui a distribuição de massa esteja 
mais concentrada no interior da estrela, aumentando da superficie para o 
centro, este efeito deverá existir, e pode afectar não só o movimento do 
plasma, mas também a estabilidade da estrela, 


A métrica das magnetosfera influencia o movimento do plasma, nao só atra- 
vés de si própria, e através do CEM dipolar, mas também através da descri- 
ção q-clássica do movimento, condições bastante diferentes do caso Min- 
kowskiano. Neste caso a descrição clássica calcula os campos pelo método 
dos potenciais de Liénhard e Wieckert (PLW), a descrição quântica da radia- 
ção assume um processo de Markov descontinuo, e o processo fecha-se 
notando-se que neste contexto se verifica o Princípio de Huyghens. Em 
espaços-G porém sabe-se não ser válido o Princípio de Huyghens e a des- 
crição clássica dos campos não pode ser feita só através dos PLW visto 
que a história anterior conta (Andrade, 1994); pelo que na descrição quân- 
tica em espaços-G a radiação de fotões já não é um efeito de Markov des- 
continuo, mas tem de reflectir a história anterior. O problema é muito mais 
complexo. 

Esta influência depende da razão R/2M.Se R/2M>2 poderemos espe- 
rar poucas diferenças da descrição em espaço-M. Contudo, se I<R2M<? 
hã efeitos fortes que podem convergir para afectar a estabilidade da estrela 
e o movimento de plasma no exterior, particularmente à superfície. Este pro- 
blema de estabilidade pode explicar porque não se observam estrelas de 
neutrões com M>2M,. 


CONSIDERAÇÕES FINAIS 


Para explicar os mecanismos de radiação de pulsares, procurámos um 
modelo da sua magnetosfera. Conseguimos separar os efeitos electromag- 
néticos dos efeitos puramente geometrodinâmicos. Vimos em detalhe o caso 
em que a magnetosfera é espaço de Minkowski, como se pode explicar a 
radiação observada em radio frequência e radiação -y e enunciámos alguns 
problemas introduzidos pela descrição geometrodinâmica. 

À investigação em pulsares tem ainda um longo caminho a percorrer. As 


condições extremas existentes em estrelas de neutrões exigem novos 
desenvolvimentos da Fisica. 
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RESUMO 


Apresentamos, aqui um método 
iterativo particularmente bem 
adaptado à resolução de siste- 
mas lineares não-simétricos, O 
qual minimiza, em cada etapa, 
a norma do vector resíduo 
sobre um espaço de Krylov. 


Palavras chave: sistemas não-simétricos, subespaços de Krylov, 


métodos descendentes, métodos de resíduo minimal. 
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1. INTRODUÇÃO 


Em Análise Numérica, um dos principais problemas computacionais é a 
resolução de sistemas lineares do tipo 


Ax=b. (1) 


Estes problemas surgem, por exemplo, quando se usa um método de diferen- 
ças finitas ou de elementos finitos, para a resolução de soluções de equações 
com derivadas parciais, na resolução numérica de problemas não lineares. 


Para sistemas lineares de baixa ordem, a aproximação clássica é feita atra- 
vés um método directo, tal como o de eliminação de Gauss. Estes métodos 
de resolução da solução de (1) baseiam-se na factorização da matriz A. 
Contudo, os sistemas lineares que surgem de problemas não académicos 
são de ordem arbitrariamente grande. Felizmente, as matrizes destes siste- 
mas têm uma estrutura especial, com muitos coeficientes nulos, pelo que 
muitas vezes os métodos directos podem ser adaptados para explorar essa 
estrutura. No entanto, em muitos casos uma aproximação directa é proibitiva 
atendendo à acumulação de erros de arredondamento e ao incremento no 
tempo de cálculo, que ela implica, pelo que a única alternativa é a utilização 
de métodos iterativos. 


2. MÉTODOS ITERATIVOS E SUBESPAÇOS DE KRYLOV 


Alguns dos métodos iterativos utilizados para resolver o sistema linear (1) 
são aplicados sobre subespaços de Krylov: com estes métodos, obtemos 
aproximações x, de A b da forma 


(2) 


onde x, e R” é uma aproximação inicial da solução de (1), r,=b-Ax, é o 
respectivo vector resíduo, e K,(r,.4) é o n-ésimo subespaço de Krylov 
gerado por r, e A, no sentido de 


KroA)=[0(A)ride Ps). 


sendo P. o espaço dos polinómios de coeficientes reais de grau inferior ou 
igual a m. Em particular, o vector resíduo obtido na n-ésima iteração pode 
ser escrito da seguinte forma 


r=b-Ax, ='P (Abr. (3) 


W7 


com 


PeP A WlO=1 (4) 


Como facilmente se depreende de (3), o objectivo dos métodos desta classe é 
em cada etapa, a determinação dum polinómio 'F, que satisfaça (3) e (4), e 
ainda tal que 1, =0. O problema inicial (1) é desta forma reduzido a um pro- 


blema de minimização, para uma dada norma ', sobre os subespaços de Krylov. 


Ir, | = min |b — À xs 
LE 4 + | 


(5) 


3. O ALGORITMO DE GMRES 


Dos sistemas lineares do tipo (1), destacam-se aqueles que resultam da dis- 
cretizaçao de operadores não-simétricos, onde por consequência a matriz A 
é não-simétrica. Para a resolução deste tipo de sistemas foram feitas gene- 
ralizações aos algoritmos dos métodos iterativos clássicos, que produzem 
aproximações aceitáveis para certas classes de matrizes não-simétricas. Em 
1986, Saad e Schultz apresentaram em [1] um algoritmo geral, particular- 
mente elegante, para a resolução destes sistemas lineares, o qual denomi- 
naram algoritmo de GMRES (Generalized Minimal RESidual algorithm). 
Sendo este algoritmo da classe referida na secção anterior, fornece-nos 
aproximações da solução do sistema (1) do tipo (2) e cuja a norma do resi- 
duo verifica (5), sem que seja necessário calcular em cada etapa x,. 


Algoritmo 1 (Algoritmo de GMRES) 
*- Amatriz A, o vector x,, o segundo membro b e a precisão £« são dados; 
e r=b-Ax, resíduo inicial; 
e om()=[hil; 
V| = | primeiro vector de base do espaço de Krylov; 
* Início de ciclo em nz: 


a 


Cálculo do n+1 vector de base 


. WD, =Av, 
. Bys=(6y ji=laaen 
. WD = - Ê. Ravi 
- hslo,| 
O 


Factorização QR da matriz A em base v. 


- Para i=l..n-l, fazer 


" Em geral utiliza-se a norma L.. 


. s=-hyr+h' 
e ho=cr-sh 
* Multiplicação do segundo membro por Q 
“ sm(n+l)=s sm(n) 
smtn)=c,sm(n) 
* Teste da norma do resíduo 
- Se |sm(n +1 )| 2€ reitere 
- Caso contrário resolva o sistema Ry=sm e prossiga 
* Fim de ciclo em n 
e Cálculo de x, 
' X, =X, + >. JM, 
* Fim do algoritmo. 


Para um sistema de dimensão N, o algoritmo (1) converge em, no máximo, 
N iterações, verificando-se o seguinte teorema de convergência (cf. [1]): 


Teorema 1 Sejam D=diag(d,.....d,,) e X uma matriz ortogonal, tais que 


A=XBX”. 


Então, temos que 


|b-Ax,|<|x||x||b- As, o": 


com 
o" = min | (max |p(d, J) 


peP,.p 


4. APLICAÇÃO DO ALGORITMO DE GMRES 


Consideremos o problema de convecção-difusão do tipo Burger's unidi- 
mensional, para exemplificar a convergência do algoritmo (1). A escolha 
deste problema não é aleatória, pois não obstante a sua simplicidade, este 
problema simula os mais importantes fenómenos computacionais da diná- 
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mica dos fluidos. Formulemos então o problema seguinte sobre o conjunto 
aberto 10, 1[: 
ox Or% 


- -v— = É em /0.l. 

arg? SM (6) 
x 0)=1, (7) 
f=0, (8) 


onde a>0,u>0e fé uma função real de variável real tomada de forma a 
que a equação (6) admita solução. 

Após discretização de (6), usando um método de diferenças finitas descentra- 
das, obtemos para x positivo, o seguinte sistema de equações não lineares: 


o e a 
é o ms + a” == E R 
As As (9) 
3 =; 
t 4 e] = 0, 


| o Ud 
Na e f =f(x,). 
Usando o método de Newton, para linearizar o sistema (9), obtemos o 
seguinte algoritmo: 


de incógnitas (x) |, com parâmetro de discretização As = 


Algoritmo 2 (Método de Newton multidimensional) 
ft. x" ee dados, 
2. para k>() resolver o sistema linear 
(= Ri) 
de MO =X'=P', 
4. iterar enquanto se verifique [rt x ) >E 
Aqui, x' e Rºº, 


W — Ri | 


Misa CE gm TATO 
REA AS : f, 


e J,(X) é a matriz Jacobiana de F em X, 
OÓbtemos desta forma uma matriz J,(X) tridiagonal não-simétrica 


+ «UR () 
"PR : 
| O a, d 
JAX)= fã 
', Ú) 
b, 
() Ud; dl, 
onde 
2ax, - 2y 
É e 
E 


p=1,.N=] 


Fazendo 


F(s)=4as' —-d4us+2v (10) 


a solução exacta do problema (6-8) será 
x(s)=1-=s”. (11) 


Usando o algoritmo (2) com 4s=0.1, a=| e v=1 obtemos uma solução 
aproximada de (11), cujo erro associado converge quadraticamente para 
zero. Na resolução do sistema do passo 2, do algoritmo 2, usamos o algo- 
ritmo de GMRES (1) o qual em cada iteração de (2) converge antes de 9 
iterações, como previsto teoricamente, com um comportamento similar ao 
exposto na Figura. 1. A eficiência deste algoritmo é superior quando utiliza- 
mos uma matriz de pré-condicionamento C (X) para a matriz J,(X), prin- 
cipalmente quando v>a. 


0.01 


10º + à 4 4 ss e 


O 1 EFA RDIO TÃEB 
NUMERO DE ITERAÇÕES 


Fig. 1 O algoritmo converge antes de 9 iterações, como teoricamente se previa. 


O algoritmo de GMRES fornece em cada etapa, e sem cálculos adicionais, 
a factorização QR da matriz A em base XK,(r,.4), utilizando rotações no 


plano. Obtém-se também a norma do resíduo da aproximação «x , sem que 
seja necessário calcular esta aproximação. Este facto constitui uma das 
grandes vantagens deste algoritmo, além de fornecer um óptimo critério de 
paragem. 


Não obstante o elevado número de vectores que o algoritmo 1 utiliza em 
memória central, a sua programação não apresenta dificuldades especiais, 
havendo somente a notar que as matrizes Q e R têm de ser armazenadas 
em memória central com especial atenção, atendendo a que são matrizes 
triangulares, construidas coluna a coluna, cuja dimensão aumenta em cada 
iteração do algoritmo. 
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RESUMO 


Neste trabalho descreve-se a 
metodologia e as técnicas 
específicas de calibração de 
analisadores automáticos de 
emissões industriais com relevo 
particular para os medidores de 
emissões de partículas totais 
em suspensão da indústria de 
fabrico de pasta de celulose, 
sendo analisadas as implica- 
ções do cumprimento da Porta- 
ria N.º 286/93 sobre a quali- 
dade do ar. 


A calibração de opacimetros 
resulta, devido aos parâmetros 
e variáveis em jogo, numa rela- 
ção multivariável, o que não se 
encontra de acordo com a cali- 
bração bidimensional que é 
efectuada na prática. 


Para o caso de uma fonte 
emissora da industria de produ- 
ção da pasta de celulose, que 
está sujeita a grandes varia- 
ções operacionais, que é a 
chaminé da caldeira de recupe- 
ração, são analisados em por- 
menor os efeitos ponderais de 
cada variável em jogo, como 
sejam, a humidade e a tempe- 
ratura do gás, desvios ao isoci- 
netismo, distribuição da granu- 
lometria das partículas e 
transmitancia caracteristica dos 
opacimetros. 


A análise multivariável de uma 
série relativamente grande de 
dados obtida ao longo dos 
anos nas unidades fabris nacio- 
nais, conduz a uma influencia 
de cerca de 98% de peso da 
transmitancia sobre as restan- 
tes variáveis, com um nível de 
significancia superior a 0,99, o 
que aponta para uma validação 
das calibrações bidimensionais 
efectuadas. 
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1. SITUAÇÃO DO PROBLEMA 


Com a entrada em vigor da Portaria N.º 286/93 para todas as instalações 
industriais existentes, o que ocorrerá em Novembro de 1995 como conse- 
quencia do disposto no Decreto-Lei N.º 352/90, coloca-se mais premente- 
mente o problema da calibração dos analisadores em contínuo o que terá 
que verificar-se obrigatoriamente uma vez por ano para as instalações 
abrangidas pelo Anexo VII da referida Portaria, ou seja, para todas as insta- 
lações que tem taxas de emissão (a taxa de emissão é igual ao produto da 
concentração do poluente pelo caudal de gás) superiores aos limites que 
constam desse anexo (1). 


De facto, para uma chaminé industrial abrangida pela Portaria N.º 286/93, 
ou seja, uma fonte estacionária fixa, podem ocorrer duas situações distintas: 
i) a concentração de poluente emitido encontra-se em conformidade com o 
valor limite aplicável (geral ou sectorial) sendo inferior a este. 


ii) a concentração de poluente emitido não se encontra em conformidade 
com a legislação sendo superior ao valor limite 


Conjugando esta informação com o caudal de gás emitido (ou seja poluente 
e gás veiculante) obtem-se então a taxa de emissão e novas situações se 
podem dar: 


a) verifica-se i) mas a taxa de emissão é superior ao especificado. 


b) verifica-se i) mas a taxa de emissão é inferior ao especificado 
c) verifica-se ii), o que implica a tomada de acções correctivas como segue: 


c 1) depois das acções correctivas já se verifica i). 
c-2) mesmo depois das acções correctivas ainda se verifica ii). 


Em termos de monitorização de acordo com os termos da legislação tem-se 
então que: para os casos a) e c-2) a legislação obriga à instalação de medi- 
dores em contínuo do poluente em causa, enquanto que para os restantes 
poluentes apenas obriga à monitorização pontual efectuada 2 vezes por ano. 


A monitorização pontual terá que ser efectuada com equipamento adequado 
e segundo os métodos USEPA ou as correspondentes normas portuguesas. 


Mesmo no caso em que existam medidores em contínuo, terá que ser efec- 
tuada, anualmente, pelo menos, uma medição pontual que servirá de cali- 
bração para os analisadores automáticos. A este respeito ver o fluxograma 
que se apresenta na figura 1. 
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Análise da fonte 
emissora: 


Cp > 
valorlimite 
? 


Acções 


Cone. poluente : Cp Correctivas 


Caudal PTN 


Cálculo da 
taxa de 
emissão : Tp 


Monitorização 
pontual ; 


em contínuo + 
calibração anual 


2? vezestano 


Fig. 1 Fluxograma indicativo do sistema de monitorização a adoptar segundo a Portaria 
N.º 286/93 


No que diz respeito à calibração dos analisadores automáticos, esta far-se-a 
sempre com base num método absoluto USEPA ou equivalente uma vez 
que estes métodos são muito mais exactos e menos sujeitos a interferencias 
num ambiente industrial agressivo do que qualquer analisador automático. 


Na verdade, enquanto o método absoluto por via húmida (USEPA ou equiva- 
lente) está geralmente bem estudado e quantificado em termos de interferen- 
cias, conduzindo a resultados representativos se forem observados os res- 
pectivos preceitos; qualquer método instrumental não é absoluto pois carece 
de calibração até mesmo porque os equipamentos de medida são, na maior 
parte dos casos, “standard” e não adaptados especificamente à medição e 
fonte em causa além de se basearem na medida de um sinal electrónico que 
é sempre falível uma vez que está sujeito a um grande número de interfe- 
rencias físicas e químicas assim como a variação do sinal face a variações 
das condições ambientais e/ou operatórias na fonte emissora. 


No caso dos analisadores de gases a calibração é feita a partir de tubos de 
permeação ou com garrafas de gases padrão, sendo, em geral a primeira 
técnica mais fiável do que a segunda uma vez que, para muitos gases a 
concentração das espécies altera-se por reacção química ou decai em con- 
sequencia da variação da pressão na garrafa. 


Encontram-se já publicadas, ou em curso de publicação, diversas normas 
ISO que referem procedimentos exaustivos (com mais de 20 ensaios com- 
parativos para cada concentração) para a calibração dos analisadores de 
gases sendo o método padrão o método de determinação do poluente por 
química húmida, segundo método USEPA ou equivalente. 


No caso dos analisadores de partículas o procedimento de calibração é 
mais complexo uma vez que, neste caso, não existem nem garrafas de 
gases padrão nem tubos de permeação aplicáveis a esta situação. É a este 
tipo de analisadores que nos referiremos em seguida. 


2. ANALISADORES AUTOMÁTICOS DE PARTÍCULAS TOTAIS EM SUS- 
PENSÃO 


Os analisadores automáticos para partículas totais em suspensão geral- 
mente utilizados na indústria são de dois grandes tipos: os analisadores B e 


ABSTRACT 


This paper describes the 
methodology and specific tech- 
niques for calibrating automatic 
on-line industrial emission 
analysers, mainly equipments 
for measuring total suspended 
dust analysers installed in pulp 
milis within the scope of Porta- 
ria N.º 286/93 on air quality. 


The calibration of opacimeters 
is a multi-parameter relations- 
hip instead of the bidimensional 
calibration which is used in 
industrial practice. For a statio- 
nary source from a pulp mill 
such as the recovery boiler 
stack, which is subjected to sig- 
nificant variations, the effects of 
parameters such as the humi- 
dity and gas temperature, devi- 
ations on isokinetism, size 
range of particles and characte- 
ristic transmitance of equipment 
are analysed. 


The multi-variable analysis of a 
considerable set of data leads 
to an estimate of about 98% of 
equipment transmitance over 
the other parameters with a 
level of significance greater 
than 0,99 which is a validation 
of the bidimensional practical 
calibrations. 


os opacímetros a que correspondem 2 princípios de funcionamento distintos. 
O primeiro tipo utiliza uma fonte radioactiva emissora B que é deflectida por 
partículas em suspensão no gás veiculante. Se bem que este tipo de anali- 
sadores conduza, em geral, a resultados mais exactos do que os opacime- 
tros, estes últimos são geralmente mais utilizados na indústria o que é 
devido ao elevado custo de investimento e instalação assim como elevados 
custos de manutenção, uma vez que requerem intervenções frequentes e, 
eventualmente, a substituição periódica da fonte radioactiva. 


Por outro lado, o opacimetro tem custos de investimento, instalação e manu- 
tenção mais baixos. Baseia-se no princípio da absorção e/ou deflecção da 
radiação luminosa de um dado comprimento de onda pelas partículas em 
suspensão no gás veiculante que atravessa o percurso dessa radiação 
(quer este seja diametral ou angular). 


Na prática os equipamentos deste tipo não medem o teor de partículas mas sim 
transmissão ou absorção de luz, o que se relaciona com a opacidade, sendo 
esta a razão pelo que estes equipamentos se designam por opacimetros. 


A extinção da luz é, de facto, uma propriedade dos aerossoís, quer seja 
medida dentro de uma chaminé ou de uma conduta, e pode ser expressa 
em termos de percentagem da luz transmitida (T = transmitancia), Op (opa- 
cidade) ou um número em escalas de densidade de fumos tais como Rin- 
gelmann, Bacharach, etc. 


A relação entre a transmitancia da luz e a opacidade da pluma de gases é 
dada por: 


Note-se que a opacidade medida no interior da fonte emissora não é, 
necessariamente, igual à opacidade da pluma de gases fora desta. As dife- 
renças entre estas duas opacidades têm a ver com as diferenças existentes 
entre as dimensões da pluma dentro e fora da fonte, diferenças nas condi- 
ções de iluminação e variações na humidade dos gases causadas pelo arre- 
fecimento da atmosfera exterior (2). 


Como se sabe, a transmissão de luz através de um aerossol pode ser des- 
crita pela lei de Lambert-Beer: 


I/1, =exp(-bL) (eq. 2) 


sendo: 


1, = fracção de luz incidente transmitida 
b = coeficiente de extinção 
= comprimento do percurso de luz 


O coeficiente de extinção b para um aerossol de diversos tamanhos desde o 
raio r1 a 12 é dado por: 


b=[O.(rLALm)a(r) n(r)dr (eg. 3) 


em que: 


O... = factor de eficiencia de extinção da luz para uma partícula de raio r, 
índice de refracção m, sendo o comprimento de onda da radiação | 

a(r) = área projectada da partícula (xr? para partículas esféricas) 

n(r) = função de distribuição dos tamanhos das partículas 


O coeficiente de extinção pode ser calculado através das equações de Mie 
(3) considerando partículas esféricas ou de qualquer outra forma. 
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A lei de Lambert-Beer pode ainda ser escrita na forma de uma relação entre 
a transmitancia ///, e a concentração mássica de partículas W, usando O 
parâmetro K definido por Pilat e Ensor (4) como segue: 


K = (volume total de partículas / volume de gás veiculante )tomn3/m3)/ 


(coeficiente de extinção da luz un —[) (eq. 4) 
ou seja: 
K=W/pb (eq. 5) 
com: 
p = densidade média das partículas (g/cm3) 
W | = concentração mássica das partículas (g/m3) 


Substituindo b, como dado pela equação 5, na equação 2, vem: 


11, =exp(LW/Kp) (eg. 6) 


e resolvendo esta equação em termos da concentração mássica W, vem: 


W=-(Kp/L) In(1/t,) (eg. 7) 


O parâmetro K pode ser calculado utilizando a teoria de Mie desde que a dis- 
tribuição de tamanhos e a forma das partículas seja conhecida. Pilat e Ensor 
(5) efectuaram cálculos de K usando distribuições de frequencia de tamanhos 
tais como log-normal e a teoria de Mie. A figura 2, cf. Carssen et al. (6) apre- 
senta os cálculos de K como função do raio médio das partículas rgw com 
vários valores do desvio padrão geométrico og para um índice de refracção 
das partículas de 1,5 e um comprimento de onda da luz de 1,2 mm. 


Este índice de refracção de 1,5 considera-se ser o típico das particulas emi- 
tidas por uma caldeira de recuperação da indústria de pasta de celulose, 
conforme (7), situação esta que iremos abordar em seguida. 


O comprimento de onda escolhido é também típico de diversos opacimetros 
industriais. Assim e, segundo Larssen (6) o valor verdadeiro de K para partí- 
culas tipicamente emitidas por caldeiras de recuperação deve encontrar-se 
entre os limites referidos na figura 2. 


Desvio padrão geometnco, TE 


Parumetro E 
comiímo 


Wc e lr 4 ara am: SME 


Raio geométrico mássco medio (Rgw), pm 


Fig. 2 Coeficiente K: volume específico de particulas/extinção. como função dos parame- 
tros de distribuição de tamanhos para um aerossol de partículas esféricas de índice de 
refracção 1,5 e um comprimento de onda de 1,2 mm, segundo Larssen (6). 


3. CALIBRAÇÃO DE OPACIMETROS 


Desde 1987 que o ISQ tem adquirido experiencia neste campo uma vez que 
procedeu já á calibração de todos os opacimetros existentes na indústria 
pesada nacional, a saber, cimenteira, siderúrgica e papeleira; com um espe- 
cial destaque para esta última, uma vez que desde a entrada em vigor do 
contrato de programa entre a SEARN e a ACEL, sobre estas unidades 
industriais recairam, desde então, as obrigações que decorrem da Portaria 
N.º 286/98. 


De facto, das várias fontes, tipos de indústria e unidades fabris em que 
estes trabalhos foram efectuados, o caso mais complexo e, como se verá, O 
das caldeiras de recuperação das unidades de produção de pasta de celu- 
lose. 


Conforme já se viu, nas equações 3 e 7, a transmitancia ou a opacidade 
são extraordinariamente dependentes da natureza e tamanho das partículas 
emitidas, o que, por sua vez, depende também das condições processuais 
de operação tais como a temperatura, velocidade, humidade e composição 
dos gases. 


Na prática, quando se faz a calibração de um opacimetro não é possivel 
tomar em linha de conta todas estas variáveis. Geralmente a calibração de 
opacimetros é efectuada de acordo com a Performance Specification 1 da 
USEPA, ou seja, admitem-se constantes todas estas variáveis e a resposta 
electrónica dos aparelhos (em mA ou %) é representada versus o teor de 
particulas em suspensão medido com um método absoluto como, por exem- 
plo, o USEPA 5 (8). No mínimo efectuam-se duas medidas correspondentes 
a duas condições de operação distintas, alterando, por exemplo, o número 
de campos activos do sistema de despoeiramento instalado a montante da 
fonte emissora. Considera-se ainda um terceiro ponto fazendo corresponder 
uma opacidade de 0% a uma concentração nula de partículas. 


Obtem-se uma representação gráfica como a que se apresenta na figura 3. 


Sinal do 
opacimetro 
(mY ou má), 
transmitancia 
ou opacidade 


io q “e 
Concentração de 
partículas, mg/Nm3 


Fig. 3 Exemplo de “recta” de calibração com 3 pontos 


Neste caso, a relação (aproximadamente) linear assim obtida pode escre- 
ver-se como: 


F=(12-r5MC2-Cl) (eg. 8) 


que se considera ser a “recta” de calibração. 
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Naturalmente, é sempre aconselhável dispor de um numero maior de pontos 
e calcular a regressão, por exemplo, pelo método dos mínimos quadráticos, 
o que aumentará consideravelmente a representatividade da calibração 
desde que o coeficiente de correlação seja aceitável (9). 


Conforme já se fez notar, este procedimento pressupõe que a opacidade é 
apenas dependente da concentração mássica de partículas mantendo-se 
constantes todas as outras variáveis atras referidas. 


Se bem que se possa admitir, sem grande erro, que isto acontece para 
determinadas fontes tais como chaminés de fornos de cimento e chaminés 
de fornos de fusão de aço em que, geralmente, as condições operatórias 
não variam muito, não variando também a natureza das matérias-primas e a 
natureza dos combustíveis, o que origina que seja pouco variável a forma e 
as dimensões das partículas emitidas, o mesmo já será mais dificil de acei- 
tar relativamente ás chaminés de fornos de cal e, principalmente, nas cha- 
minés de caldeiras de recuperação do processo kraft. 


De facto, nestas últimas fontes, os regimes de funcionamento são, como se 
sabe, muito variáveis, de acordo com as necessidades processuais (10) o 
que origina que sejam utilizados por vezes combustíveis mistos diversos 
(fuelóleo+lixívia negra com teores diferentes de sólidos+gases mal-odoro- 
sos) em situações diferentes de operação (composição percentual, tempera- 
tura, humidade e velocidade dos gases) pelo que, em rigor, a calibração 
atrás descrita não deveria ser bidimensional mas a n dimensões, tantas 
quantas as variáveis envolvidas o que seria de todo impraticável. 


Assim, e mesmo com toda a imprecisão associada a este processo, opta-se, 
geralmente, por aceitar a calibração conforme referida na Performance Spe- 
cification 1 da USEPA. 


isto origina que nem sempre se obtenha uma boa concordancia entre as lei- 
turas consecutivas dos opacimetros entre duas calibrações, uma vez que as 
condições de operação não são, de facto, reprodutiveis. 


4. INFLUÊNCIA RELATIVA DOS PARAMETROS NA CALIBRAÇÃO DE 
OPACIMETROS 


Com os resultados de determinações efectuadas em chaminés de caldeiras 
de recuperação em diversas unidades fabris, e com o objectivo de calibra- 
ção dos respectivos opacimetros, procurou-se determinar a influencia rela- 
tiva dos vários parametros atrás referidos na variabilidade da regressão 
obtida para calibração. 

As determinações de partículas totais em suspensão foram efectuadas 
segundo o método USEPA 5, e os parametros de escoamento foram medi- 
dos segundo os métodos USEPA 1 a 4, utilizando equipamentos Andersen 
segundo USEPA (12). 

Os opacimetros eram, na sua maioria Sick, de diversos modelos e encontra- 
vam-se em bom funcionamento, tendo sido objecto de manutenção e lim- 
peza antes dos ensaios. 

As amostragens pontuais foram efectuadas em locais escolhidos de 
acordo com a norma NP 2176 e durante cada amostragem efectuou-se o 
registo dos sinais lidos nos opacimetros. A determinação da granulometria 
das partículas foi efectuada por microscopia óptica segundo o método des- 
crito em (10). 

De acordo com as equações anteriormente apresentadas foi delineada a 
seguinte equação de correlação entre as diversas variáveis em jogo: 


W=AlIn(I/1,)+BT3+C fH,0+Ddgw+ E A+ F+e 


A aF = coeficientes de regressão 


Ts = temperatura do gás 

fH.O = fracção de humidade 

DI = afastamento do isocinetismo da medição em relação a 100% 
o = variação aleatória 


Para efectuar a analise de regressão multivariáveis utilizou-se um programa 
de calculo automático com base no algoritmo descrito em (11), tendo-se 
obtido os valores que constam da Tabela 1 


Coeficiente de Variação 
regressão percentual 
múltiplo (*) relativa (%) 


Nível de 
significancia 


Variável Coeficiente 
independente de regressão 


Transmitancia 
[in(l/0)] 


| Diâmetro das 
| particulas (dgw) 


Desvio do 
isocinetismo (Al) 


Temperatura do gas 


Humidade do gás 
(fH20) 


(*) calculado após a adição de cada variável 
Tabela 1 Análise de regressão multi-variável 


De facto, com base nestes valores, pode concluir-se que a transmitancia 
representa, em grande maioria, o parametro com mais peso na correlação 
(98,3%) quando comparado com os restantes 1,7% que é a soma do peso 
de todos os outros factores. 


Tem-se assim que, na gama de concentrações de partículas estudada e nas 
condições dos ensaios efectuados para as fontes em causa, a variação das 
condições operatórias (que, não foi, de facto, significativa durante os 
ensaios) é inferior a 2%. 

Nestas condições, o tipo de calibração efectuada é significativo desde que o 
opacimetro esteja equipado com um dispositivo que lhe garanta uma intensi- 
dade constante da radiação emitida, assim como da radiação recebida e/ou 
do quociente 1/lo, o que também se consegue com boa e frequente limpeza 
das lentes do aparelho, 
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ABSTRACT 


This paper presents a summary 
of the main achievements of 
the WAVEMOD project. This 
project aims at developing the 
methodology to combine data 
from different sources so as to 
reduce as much as possible the 
duration of the measurement 
campaigns that support the 
data base upon which the pro- 
babilistic models are construc- 
ted. It also aims at developing 
probabilistic models that are 
specifically applicable in coastal 
waters. The paper describes 
the field experiments aiming at 
calibrating data and describes 
the calibration of numerical 
models and remote sensed 
data. Finally the developments 
in spectral models, ARMA 
models and long-term models 
are described covering the dif- 
ferent time scales of interest in 
the description of the wave cli- 
mate, 
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1. INTRODUCTION 


The aim of the WAVEMOD project is the development of a methodology for 
obtaining probabilistic models that describe the statistics of the wave envi- 
ronment in coastal sites. The development of coastal activities requires a 
relatively accurate knowledge of the wave conditions to be expected during 
the planned life of those activities. Since predictions are required over large 
periods of time in the future, only probabilistic statements have usefulness to 
support design or operational decisions. 


Most of the probabilistic models that are available were developed to des- 
cribe the wave climate and the sea states in ocean and offshore areas which 
are characterised by deep water conditions and inexistance of complex wave 
and current fields resulting from the topography of the coast or of the bot- 
tom. In this project, consideration is given to the applicability of different 
models to the conditions available in coastal waters [1-4] (see fig. 1). 


In general, probabilistic modelling can be divided in the choice of the appro- 
priate probabilistic model and on the assessment of the value of its parame- 
ters for the data under consideration. The confirmation of the adequacy of 
the model and the quantification of its parameters is made through the 
analysis of measured data on the site under consideration. 


Fig. 1. Phenomena Generating the Coastal Waves [4] 


Since for many applications long-term models are required, it is necessary to 
build-up a relatively long series of measurements before it can be an ade- 
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quate basis to fit models for predictions and extrapolations. This project deals 
also with the problem of constructing an adequate data base by using data of 
different sources in order to reduce the length of the measurement period that 
is necessary. Thus satellite data and numerically generated data are conside- 
red in order to calibrate their predictions with measurements data. The cali- 
bration procedure is supported by two short-term measurement campaigns. 


2. FIELD EXPERIMENTS 


In order to ensure that specific problems related with different local condi- 
tions will be identified and dealt with, two sites have been chosen with diffe- 
rent bottom slopes, fetch and climatological conditions. They are in the 
Atlantic coast of Portugal and in the Northern coast of Crete in the Mediter- 
ranean [5]. 


An area on the exposed west coast of Portugal, just north of the city of 
Figueira da Foz, was identified for the Atlantic field experiment. The area 
was chosen due to its simple topography, with a relatively wide continental 
shelf, gently sloping bottom and depth contours nearly parallel to the coast. 
The wave climate was also fairly well known through long term wave measu- 
rements being made by the Portuguese Hydrographic Institute, which also 
assisted with buoy deployment. 


Fig. 2. Map of the Portuguese Field Experiment area. Location of Directional Waveriders 
are indicated by DW, WR for the Waverider, and PHI for the Portuguese Hydrographic Ins- 
titute's directional buoy. The Topex/Poseidon ground track referred to in the text is indica- 
ted by the dashed line. 


The wave climate in this area is swell dominated as a result of the frequent 
storms initiating further to the north. Local winds are predominantly northerly 
to north westerly. Wave directions are as a result predominantly north wes- 
terly although south-westerly storms do occur occasionally. 


The field experiment consisted of the deployment of 3 wave buoys along a 
line almost perpendicular to the coast from October 1993 to May 1994. 
Directional Waveriders were placed in 20m and 72m water depth and a non- 
directional Waverider was located at 50m. The Waverider was moved to 
12.5m depth in March 1994 to provide data from shallow water. The directio- 
nal buoy operated by the Portuguese Hydrographic Institute was located 
close to the 100m contour a little to the south of the buoy array and provided 
data for the latter part of the experiment. 


All the buoys transmitted directional spectra by way of the Argos satellite 
system whilst was data was stored internally in the buoys. In total, 180 days 
of directional spectra are now available from the 70m location, 150 days 
from 20m, and 90 and 80 days respectively from the Waverider at 50m and 
12.5m depth. A composite time series of significant wave height, peak wave 
period and wave direction at the spectral peak is shown in Fig 3 for the 
duration of the experiment. 
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Fig. 3. Composite time series plots showing a) signicant wave height, b) peak wave period, 
and c) mean wave direction at the peak wave period from the directional Waveriders for 
the duration of the field experiment at Figueira da Foz [5]. 


The Mediterranean site was chosen, on the north-western coast of the island 
of Crete after considering potential sites in the Hellenic seas. A small bay, 
exposed to the north and situated on the western side of the town of 
Rethymnon was identified. 


The choice was due to several factors: relatively long fetch in the Aegean 
sea, topography, sea conditions, wind conditions, and logistic facilities. The 
topography of the site chosen is not very complex and the slope is around 
2.5%. As result, diffraction and reflection are not too important. 


The site is essentially exposed to waves generated by the different local 
wind regimes with westerlies dominating. Northerlies are more common in 
summer. North easterlies also occur as a result of low pressure in the eas- 
term Mediterranean. Although summer and winter climatology are quite diffe- 
rent, severe sea conditions (Hs>4m) near shore may occur at any season. 
Wave directions at the site are predominantly north west to north east. 


Three wave buoys were deployed along a line approximately perpendicular 
to the coast. Two directional waveriders were located at 10 and 100m 
metres water depth and one non directional waverider buoy was located at 
20m depth. The campaign started on the 1st February 94 and continued until 
November 94. Both winter and summer sea conditions were covered. Three 
current-meters were also deployed near to each buoy at 8.5m below sea 
surface in order to measure the mean current velocity. The buoys transmit- 
ted raw data by means of VHF to a land station situated in a building of the 
Municipality of Rethymnon just above the site. 


3. CALIBRATION OF NUMERICAL MODELS 


The aim of this part of the project was to study the performance of hindcast 
models, as compared with measured data. The state of the art of wave 
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modelling is now in the so called models of third and even fourth generation, 
which rely heavily on the work of the WAM group [6]. These models are 
used operationally in the European Center for Medium Range Weather Fore- 
casting (ECMRWF) producing forecasts on a global scale every six hours. 
These predictions are made for grid points typically 100 km apart, which 
represent an adequate coverage for ocean areas but is clearly inadequate 
for coastal waters where the scale of the spatial variability requires informa- 
tion at much shorter distances. 


The difference between the second generation and third generation models 
is basically in the accuracy of dealing with the source term in the energy 
balance equation as well as with the non-linear wave interactions. However, 
the hindcast results are very sensitive to the input wave fields, which can 
have a much more severe impact on the uncertainty of the results than the 
differences between predictions of second and third generation models, as 
was shown recently [7,8]. 


In view of these results, the strategy of the work on this project was to con- 
centrate on the use of second generation wave models, eventually with wind 
fields produced by a local wind model with higher resolution. 


This choice resulted from considering on one hand that the objectives of this 
project were not the further development of these computationally heavy 
models for global predictions but instead interest was concentrated on 
models typically available at consulting companies dedicated to coastal engi- 
neering which will typically work with second generation models. Further- 
more the findings just mentioned, indicated that the accuracy of the predic- 
tions may be more dependent on improving the description of the relevant 
wind fields, using a local wind model with a fine grid second generation 
model, than a third generation wave model with coarsely spaced grid points. 


A wind model was developed [8,9] which improves the geostrophic wind 
fields by considering the effect of the boundary layer and calibration has 
been performed by using the predictions of wave fields based on them [8]. In 
this process one is automatically calibrating the wind and the wave model. 
The wave model adopted is based on the work of Holthuijsen and Boer [10]. 
tis a discrete spectral model with implicit coupled wave-wave interactions 
included for the region around the forecast point. Wave propagation is achie- 
ved with a tracing technique along straight rays radiating from the forecast 
point. The generation and dissipation of each spectral component is obtained 
by using the analytical solution of linear and exponential growth or exponen- 
tial decay per time step. 


The results obtained so far have shown that by appropriately choosing the 
height of the boundary layer reasonable results can be obtained [8]. 
Although satisfactory results were obtained, it is still planned to conduct fur- 
ther calibrations even with the wind fields produced by the ECMWREF. 


Another important information for coastal sites is the existing current, as a 
result of the different sources. This can be of interest on its own as well as 
in its interaction with the wave fields. A numerical model has also been 
developed so as to predict the tidal currents in the Portuguese coast [11]. 


4. CALIBRATION OF SATELLITE AND RADAR DATA 


Different types of satellite data are available and have been used in the pre- 
sent project. Altimeter data provide only the value of significant wave height 
and wind speed along the paths of the satellite. The Synthetic Aperture 
Radar (SAR) provides images from where a full directional spectrum can be 
determined by means of inversion techniques. The usage of SAR ranges 


from full resolution images for coastal studies to SAR wave mode data for 
global ocean coverage. Data is available ffom GEOSAT, ERS1 and more 
recently from Topex-Poseidon. 


At present two satellite altimeters are operating. These are the European 
Space Agency's ERS-1 Satellite launched in 1991 and US/French 
Topex/Poseidon which was launched in 1991. It is planned that ERS-1 will 
be phased out in 1995 and substituted by ERS-2 to be launched then. 


The exact repeat orbits vary from time intervals of one week to one month 
depending on the satellite and on the area under consideration. The tracks 
allow resolutions of 150 to 300Km which are acceptable for ocean areas but 
they are too coarse for coastal areas. In these areas one additional problem 
occurs when the satellite passes from land to sea, giving artificially high 
values of significant wave height. 


Referring to the Portuguese field experiment described in section 2, it should 
be noted that the directional waverider at 72m was located near to 
Topex/Poseidon ground track (fig.2). The Portuguese Hydrographic Insti- 
tute's buoy was similarly the located close to the ERS-1 ground track. 


Calibration studies on the GEOSAT significant wave height have shown a 
consistent downward bias of 13%, to be modified for low sea states. The 
ERS-1 produces significant wave heights that are significantly biased and 
needs to be corrected before use. However, the Topex/Poseidon appear to 
provide close to bias-free measurements [12]. Globally, it is generally accep- 
ted that if the global altimeter data is quality controlled and is corrected to 
remove systematic biases then accurate data is obtained. 


An example of the use of altimeter data to construct a database in the coas- 
tal area of Figueira da Foz can be found in [13]. 


Fig. 4. Inverted directional wave spectra and corresponding SAR image [5]. 


The Synthetic Aperture Radar produces images of the ocean surface with 
very high resolution which need to go through a complicated post-processing 
before extracting quantitative wave observation. The main algorithm is the 
non-linear ocean spectral inversion [14, 15]. The transformation yields an 
analytical expression corresponding to a given wave spectrum. The inversion 
procedure corresponds to finding the ocean wave spectrum such that when 
the forward ocean SAR transformation is applied to it, the results transform 
to the observed SAR spectra (fig. 4). 


Despite this procedure being routinely applied, problems exist still to be sol- 
ved. The current resolution limits the observations to waves longer than 
about 40m. Furthermore the SAR image is symmetric and thus suffers from 
180º ambiguity in direction. There is also a marked alteration in the direction 
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parallel to the radar propagation. The azimuth cut-off masks waves travelling 
parallel to the radar [12]. 


The ERS-1 provides images operationally. One image covers a square of 
100km x 100km and can be delivered as fast product (5000 x 6300 pixels) or 
high resolution (8000 x 8200 pixels). Smaller wave mode imagettes are also 
available (400 x 600 pixels). 


Recently an experimental installation of a radar has been made close to Bil- 
bao in the Northern coast of Spain. The radar is providing images that have 
some similarities with the satellite images. In fact, the method of obtaining 
images from an X-band radar is based on the effect of the capillary waves 
produced by wind, which produce interference patterns on the radar screen. 


It is possible to obtain directional wave spectra from these images and it was 
shown [16, 17] that these can be more precise than the buoy data, as con- 
cerns the directionality properties. There are however ambiguities that need 
to be resolved concerning the direction of propagation and the absolute 
value of the variance. Furthermore, without local wind there are no capillary 
waves and no image can be obtained. 


5. SPECTRAL MODELS 


5.1. Scalar Spectrum 


Over relatively short periods of time and relatively limited ocean areas, the 
elevation of the sea surface is modelled as a zero mean stationary and ergo- 
dig Gaussian process. As such, it can be completely described by a power 
spectral density function, from which the variance of the process can be cal- 
culated. 


This function includes information about the distribution of wave energy over 
the different frequencies, which is very important to several engineering 
applications. A more comprehensive description is the directional wave spec- 
tra, which provides additional information about the directional spread of 
energy as a function of wave frequency. 


At present several theoretical models have been proposed to describe wind 
waves in developing sea states and in saturated sea states. Both of these 
models rely on the Phillips model for the high frequency tail of the spectrum, 
which is dependent on the power -5 of the frequency [18]. However some 
evidence in measurements has shown that a different dependence may be 
present and theoretical works of Kitaigorodsky and Phillips have also pointed 
in that direction. This specific question has been studied in this project by 
analysing data from the two measurement sites described in section 2. 


In fact, it was found that typically the exponent will vary from -3.5 to -4.5, 
depending on the value of Hs and Tz. A clear tendency has been shown by 
the long term data from Figueira da Foz. It is believed that this variation is 
related with the state of development of the sea state, which is clearly rela- 
ted with its parameters Hs and Tz. 


Another factor predicted by the theory is the effect of finite water depth. 
While the theory predicts that in infinite depth it should be -5, an asymptotic 
value of -3.5 occurs for shallow water. This variation has in fact been obser- 
ved in the measured data in different depth as reported in [18]. 

The findings about the high frequency tail of the spectrum have been com- 
plemented with a method to predict the evolution of the stage of the spec- 
trum with decreasing water depth, taking into account the effects of refrac- 
tion and shoaling that are associated with them [19]. 
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Fig. 5. a) Deep water inverted SAR spectra (dashed), Hs=2.8m. Data from a Waverider O 
72m water depth about 80km south of the site shown with a solid line (Hs=3.3m). b) Thir- 
teen inverted SAR spectra from the outer part of the SAR scene illustrating f” decay [18]. 


5.2. Directional Spectrum 


Systematic data collection for the assessment of the directional spectrum 
has started relatively recently in different countries. While the scalar spec- 
trum is derived from time histories of the elevation of a waverider buoy, the 
directional spectrum needs information about the roll and pitch of the buoys. 
This project has been considering the different methods presently available 
to determine the directional spectrum from wave records and to derive from 
them the directionality functions [20-22]. 


Preliminary results indicate that the basic parameters of the sea state such 
as for example the significant wave height will present non-negligible diffe- 
rences when determined from spectra estimated with the Fast Fourier Trans- 
form method or with the method of maximum entropy [22]. The latter is able 
to estimate spectra from shorter time series and is able to have different per- 
formance depending on the order of the autoregressive model adopted. 
However, there is not yet a clear criterium to determine the most appropriate 
order to adopt. 


6. ARMA MODELS OF THE TIME SERIES OF SIGNIFICANT WAVE 
HEIGHT 


The time series of significant wave height represents a non-stationary process 
in time in that it is very much affected by the seasonal trends. This is well 
depicted by autoregressive models of monthly mean values as shown in [23]. 
In fact, it has been shown that the periodic nature is so clear that one term of 
a Fourier series provides an adequate model for the seasonality involved [24]. 
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Fig. 6. Observed values and estimated model from the Fourier representation [24]. 
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Modelling these seasonal effects allows the use of auto-regressive models to 
describe the deseasonalised series which in general can be described with 
models of orders ranging from 5 to about 20 although in this latter case only 
about six terms are non-zero. 


The initial experience with these models indicate that they are not appropri- 
ate to forecast the time series because they tend very quickly to the mean 
value. However, simulated series show the same properties of the initial 
series and thus can be used for several studies that are concerned with 
other derived parameters. 


One interesting application that was developed in the present project was 
the use of autoregressive models to fill the gaps of missing values in the 
time series. It was found that for relatively small gaps it was possible to 
make predictions from both sides of the gap but a finer model with orders of 
up to 65, was necessary for the purpose of fillings gaps of data [25]. 


However for long gaps of two or more weeks, which occur more often than 
desired, this method is not applicable. Instead, it was shown that one can 
use information from neighbouring measurement locations and by using 
transfer functions, to infer the missing data. 


The method, which was presented in [25], consists of fitting high order models 
to the time series in close locations and then calculating the cross correlation 
between the residuals from the two series. The structure of the cross correla- 
tion function is determined from this analysis and then it is used to determine 
the residues at one location as a function of the residues at another. 


The example shown in figure 7 shows an initial series of Figueira da Foz 
with two significant gaps at the begining and end of the record. The middle 
series showed that the method to fill gaps worked well in the second one, for 
which there was data from La Corunha, but it was not good in the first one 
which coincided with a gap in La Coruha also. 
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Fig. 7. a) Series from Figueira da Foz (1986) with the data of November removed. b) The 
series with filled gaps (in strong) based on the data from La Corunha. c) Time series from 
La Coruna (1986) [25]. 


The examples considered have shown that this method works whenever 
both locations are subjected to the same wave climate. In fact, using two 
stations located at approximately the same distance from a reference buoy, it 
was shown that in one case a good correlation was obtained while no corre- 
lation existed for the other where local disturbances existed. 


Another approach to the modelling of time series can be found in [26] 
although applied to hindcast data from the North Atlantic. 


7. LONG-TERM PROBABILITY MODELS 


The description of the variability of the wave climate over periods of many 
years is of interest to design formulations of different types of fixed and floa- 
ting structures. The current procedure is to organise a data set with several 
years of measurements, to fit a distribution and afterwards to extrapolate it to 
very low levels of probability. 


This procedure has started to be critically examined in [27] where the con- 
cept of the underlying distribution was discussed pointing out that in general 
the current procedure would be indeed pooling together data that was from 
populations with different characteristics. This work was pursued in [28] 
where it was shown that the basic assumption that the data from successive 
years could be considered as a random sample from the same population 
was not supported in general by statistical tests of fit. 


Another aspect of interest is related with the parametric models that are 
used to describe the significant wave height data. Frequently used models 
are the log-normal distribution, when concern is the lower to moderate range 
and the Weibull distribution when a better model of high sea states is aimed 
at. However, a combination of both have also been proposed [29], as well as 
the Gamma distribution [30]. 


The choice of the probabilistic models is often justified by a the visual appea- 
rance in a plot of probability paper. Furthermore, often the fit of the distribu- 
tion i.e., the assessment of its parameters is made by the method of least 
squares which is a formalization of the method based on probability paper. 


It was shown in [27] that when using other methods to fit distributions, such 
as the method of moments and the method of maximum likelihood, different 
values are obtained for the parameters. As expected, the difference on the 
parametric values are larger whenever the quality of the fit is not good. 


The tests of goodness of fit can indicate whether a distribution should be rejec- 
ted. However, often these tests are not robust enough to select only one distri- 
bution. This was shown in [31], where different distributions based on Beta 
functions were used and not rejected by the used tests of fit. It was argued 
that whenever extrapolations were made with Beta distributions with different 
tail behaviour they provided bounds to the extrapolated values to be expected. 


The different long-term models of marginal distributions of significant wave 
height are adequate for several applications but not for all. In some situati- 
ons it is necessary to model the bivariate distribution of significant wave 
height and mean or peak period. 


The method that has been commontly used is to combine the marginal distri- 
bution of significant wave height with several distributions of period conditio- 
nal on Hs [29]. However, several difficulties occur for high sea states, in 
which few observations are available and the fit of the conditional distributi- 
ons is difficult. 


Alternatives to this procedure is the use of the Placket model, which was 
proposed recently and applied to North Atlantic data [32], or the use of trans- 
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Fig. 8. Three dimensional plots of the Kernel density estimates of the (Hs, Tz) data of 
Figueira da Foz (1986/87) and the fitted density function [33]. 


formation of variables which was developed in the present project and 
applied to data of Figueira da Foz [33]. 


8. CONCLUSIONS 


An overview has been presented of the work that has been performed roughly 
in the first part of the WAVEMOD project covering the field experiments to col- 
lect the data for calibration, the work to calibrate numerical models and remote 
sensed data and finally the new achievements as concerns the short term 
spectral modelling of sea states, the intermediate scale ARMA models of the 
sea state time series and the long term models used eventually for extrapola- 
tion to severe conditions of low probability of occurrence. 
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ABSTRACT 


It is hard to pick up a technical 
journal today that does not 
include some information on 
new electrochemical sensors. 
And even in the popular litera- 
ture we have seen articles on 
topics such as chemically sen- 
sitive field effect transistors 
(CHEMFETS) for in vivo moni- 
toring. The present article inclu- 
des a general description of ion 
selective microelectrodes 
(ISMs), being the first of a 
series of papers on advanced 
chemical devices for clinical 
applications, 
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1. INTRODUCTION 


The importance of extracellular and intracellular recording with ion selective 
microelectrodes (ISMs) is growing rapidly in several areas of experimental 
physiology [1-4]. These are particularly useful in renal, gastrointestinal, and 
cardiac electrophysiology, as well as in studies of the microcirculation, and 
of membrane transport in a variety of cell types. Today the term “microelec- 
trode” implies electrodes with tip diameters that do not exceed 1-2 um, and 
very often tip diameters are even smaller than 1 um. Our particular interest 
in ion selective microelectrodes arose because of our desire to make mea- 
surements of very small ion concentrations of chemical species with rele- 
vance in clinical chemistry. ISMs have been combined with ion selective 
field effect transistors (ISFETs), and constitute a first step to the develop- 
ment of more advanced devices for continuous in-vivo monitoring of certain 
electrolyte species. 


2. KINDS OF ELECTRODES 


There are three kinds of microelectrodes that have been used so far in biolo- 
gical studies and these are solid state, glass membrane, and liquid mem- 
brane microelectrodes. At present, the latter type is the most widely used. 
Liquid membranes can be based on ion-exchangers [5, 6] or on neutral carri- 
ers [7-10]. 

Figure 1 illustrates a liquid membrane ISM. It is composed of a glass microe- 
lectrode, the tip of which is filled by an ion selective liquid membrane. The 
contact to the back of the liquid membrane is achieved by means of a filling 
solution which contains the ion for which the membrane is selective, as well 
as a species which defines the potential of the internal contact with respect 
to the filling solution. Generally, the internal contact is achieved by means of 
an AgCl-coated Ag wire, the potential of which is determined by the pre- 
sence of Cl ions in the filling solution. 

lon selective liquid membranes are hydrophobic organic liquids; on the other 
hand, a clean glass surface is hydrophilic, which means that the liquid mem- 
brane will be displaced from the microelectrode tip by water if the glass sur- 
face is not treated to become hydrophobic [11]. The major component of 
glass being silicon dioxide, the oxygen atoms at the surface can be present 
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Fig. 1. Schematic representation of a liquid membrane ion selective microelectrode. 


in one of two forms as shown in Figure 2. Siloxane oxygens are bonded to 
two silicon atoms while a silanol oxygen is bound to only one silicon atom. 
The proton of the silanol can dissociate and leave a negatively charged oxy- 
gen on the surface. These negative sites can then interact with cations con- 
tained in the solution with which the glass is in contact and the glass surface 
can behave as a cation electrode. It is therefore important to mask as many 
of the silanol groups as possible in order to optimize electrode function. This 
is done by means of a silanization treatment which can be conducted either 
in liquid or vapour phase. Figure 2 shows an ideal situation in which the sila- 
nol groups are inactivated by reacting with trimethylchliorosilane. 
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Fig. 2. Conversion of a siloxane group to two silanol groups by hydration, and reaction of 


the silanol groups with trimethylchiorosilane. 


3. MEASURING CELLS 


A typical measuring cell is shown in Figure 3. By selectively transferring the 
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Fig. 3. Schematic diagram of ion selective microelectrode measuring circuit and cell 
assembly. 


measured from the sample solution to the liquid membrane phase, a potential 
difference is generated between the internal filling solution and the sample 
solution [12]. Ideally, the potential E of the cell is given by the Nernst equation: 


E = E, + sloga, (1) 
However, because of the presence of interfering ions, the potential of the cell 


is described by the successful semi-empirical relationship developed by 
Nicolsky and Eisenman [13, 14]: 


; = E, +slog . + RS (a) | 


(2) 


with 
A 
- 2.303 RT o 59.16 mV (at 25º €) 
Z, H Z, (3) 
and 
E, = E, + E, + Es (4) 
where 


E = potential of the cell 


E. = constant potential difference including the boundary potential diffe- 
rence between the internal filling solution and the membrane. 


E, = constant potential difference consisting of the potential difference 
between the metallic lead of the inner reference electrode of the ISE 
and its internal filling solution, and of the potential difference bet- 
ween the reference electrolyte and the metallic lead of the reference 
electrode. E, is independent of changes in sample composition. 


E, = liquid junction potential difference generated between reference 
electrolyte and sample solution. E, is variable and sample solution 
dependent. 


z. a, = Charge number and sample solution activity of the primary ion i. 
Z,. a, = charge number and sample solution activity of the interfering ion j. 
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kK”” = potentiometric selectivity coefficient; K/” is a measure of the prefe- 
rence by the sensor of the interfering ion j relative to the primary 
ion 1. 
R = molar gas constant. 


+ 
H 


absolute temperature. 
F = Faraday constant. 


4. ELECTRICAL NOISE 


The high resistivity of ion selective liquid membranes combined with the very 
small tip diameter of ISMs results in high electrical resistance of these devi- 
ces (10º 0 - 10" Q) [1, 5, 7, 10]. This can cause difficulties with electrical 
noise. 


Even if it is not generated by the system, electrical noise can be present due 
to capacitive, inductive, and electromagnetic coupling of the system with an 
external noise source [15]. In the case of an ISM, the most bothersome 
source of noise is due to motion of an operator around the physiological pre- 
paration; in this case, the noise is due to capacitive coupling with the voltage 
due to charge build-up on the operator's body. Figure 4a shows the situation 
schematically. The reference electrode is grounded, and the ISM is connec- 
ted to an electrometer by means of a wire; C, represents the stray capaci- 
tance between the wire and ground and C., represents the 


a 
TO ELECTROMETER 


REF ISM 


Cc 


Fig. 4. a) Capacitive coupling of noise source to ISM. 
b) Equivalent electrical circuit. 


coupling capacitance between the ISM-wire assembly and the noise source 
Vc: V, |s the noise voltage developed across the ISM. Assuming there is 
no capacitive coupling with the reference electrode, the impedance of the 
reference electrode is negligible compared to that of the ISM, and the ISM is 
purely resistive, the electrical equivalent circuit for the arrangement is the 
one shown in Figure 4b. Considering this circuit, the amplitude v, of the 
noise voltage produced across the ISM is given by: 


NE = RC, V 


N I 5 


+ q [(c, + Cs) RO (5) 


This equation does not show explicitly how the noise voltage depends on the 
various parameters, however, two limiting cases can be considered. First, we 
consider the case when R is much lower than the impedance of CC. + C,: 


R <<-—-——————— 
m(C. + €,) 


In this case, Equation (5) reduces to: 


Vu = RC. V. (6) 


This clearly shows that, in this case, V, is proportional to the angular fre- 
quency «, the resistor R, the coupling capacitor C,.., and the amplitude of 
the disturbing voltage V.. A second limiting case is expressed by: 


| 


R >> ———————— 
m(C. + C) 


In that case, Equation (5) becomes: 


Vas ge di 
c TC (7) 


This is, of course, the limiting case for high frequencies. A plot of the noise 
voltage V, as a function of the angular frequency « is given in Figure 5. 


YN fo Vy =WRCçVo 
/ 


Fig. 5. Frequency response of capacitively coupled noise voltage. 


Since in the case of microelectrodes, one has no control on the noise sour- 
ces and their characteristics, or on the magnitude of the resistance R, the 
only options available for the reduction of the noise voltage V, are shielding 
or reduction of the length of the connecting leads, which will result in reduc- 
tion of the coupling capacitor C.. The latter can be achieved by directly 
mounting the ISM tip on top of the gate of an ISFET as it would be discus- 
sed in a future publication. 


5. FINAL REMARKS 


Intracellular ISMs using ion selective liquid membranes have a very high 
resistance because of their very small tip diameter (0.5-1.0 um). The value 
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of the resistance ranges from 10'0 1 10"Q [1], depending on the tip size 
and the type of ion selective liquid membrane. This imposes the requirement 
of a very high input impedance and a very small input current on the measu- 
ring electrometer. Any shunt capacitance across the ISM such as shielding 
of the connecting lead will increase the response time of the probe. On the 
other hand, shielding is necessary in order to reduce the influence of elec- 
trostatic and electromagnetic noise. It is generally most desirable to bring the 
electrometer as close as possible to the microelectrode. One possibility is a 
microprobe electrometer in which the preamplifier is actually part of the 
ISM's holder. A further step is to place the microtip directly on the gate of an 
ISFET acting as the amplifier, leading to the micro ion selective field effect 
transistor (MISFET). Using appropriate liquid ion exchangers (LIES), MIS- 
FETs can be developed and characterized in terms of sensitivity, selectivity, 
limit of detection, electrical response time and chemical response time. In 
order to measure chemical response time, special purpose flow injection 
analysis (FIA) systems can be envisaged. Transfer function measurements 
can also be conducted on MISFETs and interpreted in terms of the impe- 
dance of the microelectrode placed on their gate. Most of these studies are 
being of our concern for some time. This short paper provides basic informa- 
tion about the ion selective microelectrode, as a pioneer concept in this 
broad domain of chemically sensitive solid state devices for medical and bio- 
logical chemistry. 
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RESUMO 


Este artigo propõe descrever a 
evolução morfológica da cidade 
de Lisboa desde a época medi- 
eval até à actualidade. O objec- 
tivo é compreender o desenvol- 
vimento global da cidade sob o 
ponto de vista da sua organiza- 
ção espacial relacionando-o 
com as transformações ocorri- 
das ao longo da sua história, 
detectando as situações de 
continuidade e ruptura ocorrfi- 
das no seu tecido urbano e 
desenhadas no tempo. A abor- 
dagem seguida baseia-se na 
descrição morfológica dos fac- 
tos urbanos e utiliza como 
metodologia a teoria de Space 
Syntax (Hillier e Hanson, 1984). 
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1. INTRODUÇÃO 


Este artigo tem por finalidade descrever a evolução morfológica da cidade 
de Lisboa desde a época medieval até à actualidade. O principal objectivo é 
entender o desenvolvimento global da cidade sob o ponto de vista da sua 
organização espacial relacionando-o com as transformações ocorridas ao 
longo da sua história. Pretende-se detectar as situações de continuidade e 
ruptura ocorridas no seu tecido urbano e desenhadas no tempo. A descrição 
morfológica dos factos urbanos utiliza como metodologia a teoria de Space 
Syntax (Hillier e Hanson, 1984). 


O estudo é iniciado com a descrição da cidade medieval. Este primeiro 
momento de análise baseia-se em elementos fornecidos pelo mapa de 
J.N.Tinoco (1650), considerada a fonte cartográfica mais antiga e fidedigna 
da cidade de Lisboa. Este mapa evidencia a relação da cidade — compreen- 
dida entre o Bairro Alto e o Campo de Sta. Clara —- com a zona ribeirinha 
revelando ainda o tempo dos descobrimentos portugueses. Segue-se a aná- 
lise do novo traçado da cidade Pombalina (terceiro quartel do Séc.XVIII) que 
introduz uma ruptura morfológica: a cidade volta 'as costas para o rio' e pro- 
move a ligação entre as partes oriental e a ocidental, mantendo no entanto 
algumas permanências — Rossio, Praça do Comércio, Castelo, Alfama e 
Bairro Alto. A rectificação do traçado configura uma nova cidade que cada 
vez se afastará mais do rio, sua anterior principal razão de ser, passando-se 
da Lisboa do apogeu das Descobertas para a Lisboa que contém o germe 
da futura metrópole. Após a intervenção pombalina, a manutenção do sen- 
tido dos traçados urbanos confirma a interiorização de Lisboa pelo desenvol- 
vimento de linhas de fixação ao longo de dois vales: do VALVERDE (Ave- 
nida de Liberdade — Avenida da República) e do vale dos Anjos (Avenida 
Almirante Reis — Avenida Gago Coutinho). Este desenvolvimento esboça-se 
na Lisboa Romântica estruturando-se com o traçado da Lisboa de Ressano 
Garcia para posteriormente se configurar na Lisboa do Estado Novo ante- 
cessora da Lisboa actual. 


2. MÉTODO DE ANÁLISE 


A abordagem utilizada neste trabalho consiste numa aproximação racional à 
descrição do espaço, entendido como entidade morfológica fundamental no 
campo da arquitectura e do desenho urbano, inicialmente desenvolvida por 
Hillier e Hanson (1984). Traduz-se num conjunto de técnicas descritivas — 
representação, quantificação e interpretação — da estrutura configuracional 
de sistemas espaciais. 
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A estrutura configuracional é o resultado das relações estabelecidas entre 
os espaços de um sistema tendo em atenção todos os outros que o consti- 
tuem ou mais precisamente, como os espaços se articulam e relacionam 
entre si enquanto elementos de um conjunto mais lato. O sistema espacial 
constitui uma estrutura de natureza contínua definida pela articulação entre 
os elementos construídos e pelas relações de permeabilidade e visibilidade 
estabelecidas. 


Este sistema é representado gráficamente pelo mapa axial — transformação 
bidimensional feita a partir da planta — constituída pelo menor conjunto das 
maiores linhas que cobrem todo o sistema espacial e que correspondem à 
imagem de continuidade física e visual experimentada por aqueles que se 
movem ou permanecem estáticos no sistema. O mapa axial permite retirar 
informação alfanumérica de modo a quantificar as propriedades sintácticas 
da configuração em função dos níveis de acessibilidade presentes. Para 
cada linha axial são atribuídos valores numéricos que expressam a sua 
importância em relação às linhas vizinhas e a complexidade dos percursos a 
partir de cada linha em relação às restantes. Para tal, o mapa axial é trans- 
formado num grafo — o grafo axial, As linhas no mapa axial correspondem 
aos vértices desse grafo sendo que dois vértices são adjacentes se e só se 
as linhas correspondentes se intersectarem. O grafo axial contém informa- 
ção acerca das conexões de cada linha com todas as outras que constituem 
o sistema. O conceito de acessibilidade é entendido não em termos de dis- 
tância métrica mas em função do número de mudanças de direcção — pro- 
fundidade do grafo — que é necessário efectuar para ir de um espaço a 
outro no sistema. 

A medida sintáctica de maior relevância é a Integração porque define o grau 
de acessibilidade de cada linha axial em relação a toda a área de estudo 
ou, por outras palavras, porque define o modo como o sistema de percursos 
a partir de uma linha em relação às restantes é, relativamente, simples (inte- 
grador) ou tortuoso (segregador). O grau de integração ou segregação que 
uma linha axial (espaço) apresenta em relação às restantes linhas do sis- 
tema está dependente do número de mudanças de direcção necessárias, ou 
seja da profundidade da linha. Estes valores são medidos através da Assi- 
metria Relativa (RA). Esta medida exprime a relativa profundidade do 
espaço a todos os outros do sistema. O valor de RA para a linha axial i é 
calculado pela seguinte expressão: 


RA, = 2(DMA)(K-2) 


em que: 


DM, = a distância média da linha axial | em relação as restantes linhas do 
sistema, 
K = número total de linhas axiais presentes no sistema. 


Os valores de RA variam entre O e 1. Os valores mais baixos indicam que o 
espaço tende a integrar o sistema enquanto que a valores mais elevados 
correspondem a linhas axiais ou a espaços urbanos com tendência a serem 
segregados pelo sistema. 


A Assimetria Relativa é uma função do número de total K de espaços ou 
linhas axias presentes num sistema. Para efeitos comparativos entre mapas 
axiais de diferentes dimensões, os valores de RA são padronizados pela 
Assimetria Relativa da raiz de um grafo em forma de diamante com igual 
número K de espaços. Ao inverso desta medida dá-se o nome de Integra- 
ção. Hillier e Hanson (1984) tabelam RA da raíz desta forma em diamante 
em função do número K de espaços e Kruger (1990) apresenta a expressão 
matemática que traduz numéricamente estes valores. Como veremos os 


valores da integração média das linhas axiais dos mapas em análise são 
função do número K de espaços pelo que, para atenuar de forma significa- 
tiva este efeito, descreveremos a variação da integração pelo seu coefici- 
ente de variação dado pela razão entre o desvio padrão e a média da inte- 
gração de cada mapa axial. 


A Integração define uma propriedade global do sistema: O modo como os 
espaços se relacionam com os outros imediatamente adjacentes é expresso 
através de medidas locais: conectividade, controle e raio 3. A conectividade 
exprime o grau de conecção de um espaço em relação aos espaços vizi- 
nhos i.e., indica o número de espaços ou linhas axiais imediatamente adja- 
centes. O controle traduz a quantidade proporcional de acessos que um 
espaço apresenta em relação aos espaços adjacentes i.e. indica o grau de 
importância que um espaço assume em relação aos espaços vizinhos 
enquanto ponto obrigatório de passagem. Os valores mais altos correspon- 
dem aos espaços com maior controle nas relações locais. Os valores mais 
baixos correspondem a espaços com um fraco controle local. O raio 3 mede 
a integração local dos espaços i.e. o número de linhas e mudanças de 
direcção que são necessárias para ir de um espaço a todos os outros que 
se encontram a uma profundidade igual ou menor que 3. O seu valor é cal- 
culado através de fórmula idêntica ao valor da integração mas sómente são 
consideradas distâncias iguais ou inferiores a três unidades. 


Se considerarmos que A, representa a matriz de adjacência do grafo axial 
em que A, = 1 se o vértice i for adjacente ao j e, caso contrário, A,= O então 
a Conectividade C' da linha axial i será dada pela seguinte expressão: 


C', ns 2, A, 
e o Controle Ci da mesma linha i sobre as restantes será dado por: 
C:, = 2; ( A, 3 A, 


A correlação estatística entre a medida global e a local estima a inteligibili- 
dade sintáctica do sistema. Esta medida relaciona o grau de conexões de 
uma linha — propriedade local — com a posição que esta assume na globali- 
dade do sistema — propriedade global. A inteligibilidade sintáctica expressa 
o modo como o sistema espacial em questão pode ser reconhecido no seu 
todo através da leitura das suas partes constituintes. 


A compreensão da estrutura sintáctica e respectivos padrões espacias pre- 
sentes numa configuração é facilitada pela leitura dos cascos obtidos em 
função de cada propriedade. O casco é a amostra constituida pelos 5%, 
10% ou 20% (dependente da dimensão do sistema) espaços ou linhas axi- 
ais que obtiveram maior valor, numerados por ordem decrescente. 


3. BREVE ENQUADRAMENTO HISTÓRICO 


A fim de enquadrar históricamente a evolução urbana de Lisboa descreve- 
se, brevemente, a expansão e a relação com o rio contrapondo-se a cidade 
ribeirinha à cidade do interior. 


Nos finais do século XV e sobretudo no início de Quinhentos, tempo que 
corresponde afinal ao período de expansão além-mar, Lisboa cresce consi- 
deravelmente para abrigar os novos habitantes que os negócios, a constru- 
ção naval, a “carreira da Índia"(Moita 1983) atraíam. Os 60 mil habitantes do 
século XV passam a 80 mil no século seguinte para chegarem a atingir os 
150 mil em seiscentos (Vieira da Silva cit. por França 1987), números reve- 
ladores de um acréscimo considerável de dinâmica urbana. 
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As estruturas medievais são abandonadas. O rei já não se abriga atrás de 
inúteis muralhas. Constrói um novo palácio real da Ribeira ao lado dos esta- 
leiros navais, junto ao rio, determinando o futuro desenvolvimento da cidade 
que definitivamente encontrava no mar a sua razão de existir. E precisa- 
mente a alteração das funções da cidade e a localização das novas activi- 
dades que vão reger as transformações urbanas (Marques, 1994). Desen- 
volvimento ribeirinho e crescimento para o interior vão definir os parâmetros 
do desenvolvimento de Lisboa. 


O lugar de Lisboa nasce da relação terra-rio/mar. A enorme bacia que o 
Tejo forma frente á cidade é o elo desta ligação e permitiu ao português de 
quinhentos o usofruto do melhor porto da época. Já em meados de quatro- 
centos quando o centro urbano e político continuava no castelo, o foco 
comercial já se desviara em relação á baixa e eram abertas portas na mura- 
lha para a comunicação com o Tejo. Mas é com D.Manuel | que são feitas 
as grandes obras de aterros na Ribeira, incluindo a criação do Terreiro do 
Paço, que definitivamente, marcará a cidade como uma presença determi- 
nante no seu desenho — retomado na cidade pombalina — e junto ao qual 
será erguido o novo Paço Real para onde o se transfere em 1505, consa- 
grando o terreiro do paço como centro político e administrativo. A cidade 
medieval deu lugar a uma cidade moderna e comercial, aberta ao rio, ao 
Atlântico. 


A estreita relação com o rio vai manter-se até que os planos delineados no 
pós terremoto, quer para a nova Baixa de traçado iluminista que se ergue 
sobre os escombros da velha cidade, como para os estudos utópicos de 
expansão para o interior de Lisboa ( Moita 1982, p.87), anunciem clara- 
mente a intenção dos arquitectos Pombalinos da inevitabilidade do cresci- 
mento norte da cidade. 


Desenhado no remate da rigidez do traçado da Baixa, logo após a nova 
Praça do Rossio, o Passeio Público, para além de constituir um inédito pro- 
grama de equipamento verde acrescentado à nova cidade iluminista, situ- 
ava-se cronológica e topograficamente como um novo jardim localizado na 
continuação da cidade pombalina, apontando conscientemente no sentido 
do desenvolvimento urbano (França, 1987, p.138). 


Com efeito, quando em 1879 é aberta a Avenida da Liberdade sobre o Pas- 
seio Público, o boulevard rasgado para norte, cuja ideia fora lançada em 
1859, apenas retoma uma intenção dos arquitectos pombalinos que, sem 
porem em causa o corpo ribeirinho da cidade intuiram a funcionalidade de 
um eixo ortogonal cortando o Valverde e Vale Pereiro em direcção à circun- 
valação (Henriques da Silva 1989, p.21), finalmente definida em 1885. 


A Avenida da Liberdade transforma-se na “Avenida” duma cidade assim 
libertada da sua configuração tradicional, que a reconstrução pombalina 
ainda adoptara, definindo-se como eixo de penetração da cidade que se 
construía ao longo do século XX, iniciada a partir dos bairros Camões e 
Barata Salgueiro que a ladeavam, prolongada pelas avenidas Novas até ao 
Campo Grande. Paralelamente um eixo mais pequeno-burguês o da Rainha 
Dona Amélia depois Almirante Reis, prolongava a rua da Palma, segundo 
eixo de penetração ao longo de um vale que a cidade do estado Novo, a 
partir de finais dos anos 20, completaria e condensaria numa estrutura mais 
vasta, consagrada no Plano Director de De Groer (1938-1948). A Cidade do 
Estado Novo desenvolve-se para Alvalade consagrando definitivamente a 
expansão a oriente, para Belém-Restelo estruturando a zona oriental e pon- 
tuando a cidade com importantes equipamentos (IST-Alameda, Monsanto, 
Aeroporto...) 


Por isso, alguns autores consideram que Lisboa, de cidade organizada ini- 
cialmente em função do rio, enceta a partir do 3.º quartel do século XIX uma 
nova expansão, desta vez em direcção ao planalto, com uma estrutura radio- 


cêntrica ramificada a que corresponde num primeiro momento a construção 
da Avenida da Liberdade seguida das Avenidas Novas, e num segundo 
momento, considerado do Estado Novo, fechando a expansão oriental e oci- 
dental (Nunes da Silva 1987, p.48) como que definitivamente marcando o ter- 
ritório urbano da futura metrópole, isto é, o centro da área metropolitana de 
Lisboa esboçada ao longo da cidade dos anos 60 até à actualidade. 


4. DESCRIÇÃO DOS MAPAS UTILIZADOS NA ANÁLISE 


As fontes cartográficas utilizadas obedecem a dois critérios de selecção: 
fidedignidade e representavidade balizando-se os momentos mais importan- 
tes da evolução da cidade. 


A Planta Topográfica de Lisboa em 1650 por José Luís Tinoco (Fig. 1) é 
a primeira planta da cidade que, apesar de levantada em seiscentos, apre- 
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Fig. 1 Planta de Tinoco 
Planta topográfica de Lisboa em 1650 por João Nunes Tinoco 


senta a estrutura medieval. Inclui a ocupação ribeirinha, mostrando a trans- 
formação e ocupação das muralhas nesta zona. É limitada a ocidente por 
um troço do Bairro Alto e a oriente pela Graça (Cerca dos Frades de 
S.Vicente, Campo de St.” Clara, St.? Engrácia) e pela fundição de Artilharia 
junto ao rio. Permite reconhecer a Cerca Moura e a Cerca Fernandina (assi- 
naladas a vermelho), a localização das diversas ruas por freguesias e dos 
equipamentos principais. É clara a presença do Rosssio e do Terreiro do 
Paço e a sua forte ligação através da estrutura urbana orgânica caracteri- 
zada por um emaranhado de ruas medievais, acentuando o contraste entre 
Alfama, Mouraria, mesmo a Baixa, com o loteamento quinhentista, reticulado 
e fora de muralhas, do Bairro Alto. 


A Planta Pombalina (Fig. 2) reproduz a cidade do pós-terremoto, colorida, 
assinalando com banho de vermelho o que se conserva antigo e amarelo o 
“Projecto do novo Plano” delineado pelos engenheiros militares. À cidade 
barroca de D. João V, com o aqueduto abastecedor, sobrepõe-se a nova 
cidade “iluminista”. Os limites urbanos alargam-se: a ocidente até Alcântara, 
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a oriente às Comendadeiras de Santos, a norte desenha-se a penetração 
planeada do Passeio Público (primeiro equipamento verde público) acompa- 
nhada pela rua de S.José delineada até S. Sebastião, constituindo uma 
das saídas norte da cidade tal como o regueirão dos Anjos. O desenho da 
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Fig. 2 Planta Pombalina 
Planta topográfica da cidade de Lisboa na segunda metade do século XVII 


cidade é alterado com o traçado iluminado da Baixa levantada sobre os 
escombros do terremoto, usando-se uma malha reticulada que se adapta a 
outras intervenções pombalinas como o Bairro da Fábrica das Sedas as 
Amoreiras ou o conjunto da Lapa. A sul a margem ribeirinha é ampliada 
com o monumental prolongamento do T. Paço, a nova Praça do Comércio, 
que mantém uma forte ligação ao Rossio. 


A Planta Romântica (Fig. 3) representa uma planta com vista de Lisboa em 
1833 de origem inglesa. Os limites da cidade mantêm-se, alargando-se ape- 
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Fig. 3 Planta Romântica 
Planta e vista de Lisboa em 1844 


nas a ocidente a ocupação em Alcântara e a oriente até Stº Apolónia. Os eixos 
de penetração norte clarificam-se, sobretudo rua do Salitre e a rua de S.José 
contornando o passeio Público, e a sul inicia-se o aterro da Boavista. Trata-se 
da consolidação do tecido delineado no tempo de Pombal e da definição de 
uma linha de defesa circundando a área urbana desde Alcântara a Santos, 
passando por Campo de Ourique, Picoas, Arco Cego, Alto de S. João. 


A Planta de Ressano Garcia (Fig. 4) consiste na planta topográfica publi- 
cada em 1871 tendo sobrepostas as alterações feitas até 1911. Os limites 
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Fig. 4 Planta de Ressano Garcia 
Planta topográfica de Lisboa publicada em 1871, tendo sobrepostas a tinta encarnada as 
alterações feitas até 1911. 


mantêm-se em Alcântara e alargam-se até à Madre de Deus. A sul o aterro 
da Boavista e, perfeitamente definida a Estrada de Circunvalação articulada 
com a expansão norte da cidade desencadeada a partir da abertura da Ave- 
nida da Liberdade (sobre a destruição do Passeio Público). Definem-se os 
bairros Camões, Barata Salgueiro, o Parque da Liberdade e as Avenidas 
Novas que se estendem até ao Campo Grande. A ocidente, Campo de Ouri- 
que, o bairro ocidental ao Parque da Liberdade e a oriente a abertura da 
Rainha Dona Amélia, actual Av. Almirante Reis, no esteio do regueirão dos 
Anjos. Novos equipamentos pontuam a cidade (elevadores, eléctricos, gares, 
mercados, grandes armazéns) e a localização de industrias é acompanhada, 
por vezes, pela construção de vilas operárias. 


A Planta do Estado Novo (Fig. 5) define-se com as alterações efectuadas até 
1948 sobre os levantamentos de 1899-1904-1911. A área urbana corres- 
ponde exactamente à actual, isto é, dos Olivais a Algés. Inclui os aterros 
efectuados em Belém (1940) e as Gares Marítimas, os Bairros Económicos 
do Arco do Cego, Encarnação, Madre de Deus, Caselas, Salazar e ainda o 
início do Bairro de Alvalade (assinalando-se as infraestruturas viárias bási- 
cas). A Avenida Almirante Reis é concluída com a Avenida do Aeroporto (o 
novo equipamento no limite norte da cidade) passando pela Praça do Are- 
eiro, urbanizando-se toda a área da Alameda definida pelo IST em articula- 
ção com as Avenidas Novas, assim se ligando à cidade esta nova zona ofi- 
ental de expansão. À ocidente o plano do Restelo articula-se com o novo 
equipamento verde, o Parque de Monsanto. Paralelamente ocorrem interven- 
ções mais pontuais que se situam em linhas axiais bem definidas: Avenida 
Oriental ao Parque, D. João V, Paiva Couceiro, etc., imprimindo à cidade 


56 


alguma disciplina urbana que se liga à abertura de grandes eixos viários arti- 
culados com as saídas da cidade por auto-estradas (Norte, Estoril e Sintra). 
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Fig. 5 Planta do Estado Novo 
Organizada em 1899 assinala as transformações realizadas entre 1911 e 1948 


A Planta Actual (Fig. 6) baseia-se num levantamento efectuado pelos Servi- 
ços Cartográficos do Exército. Os limites da área urbana permanecem inal- 
terados relativamente à anterior. Inclui os bairros de promoção pública com 
Alvalade e Olivais concluídos, Chelas e Telheiras, e toda a expansão de 
carácter especulativo e desordenado ao longo dos eixos periféricos de Ben- 
fica, Luz, e Lumiar. Em termos de novos eixos viários assinala-se a 2.º cir- 
cular e a Ponte sobre o Tejo estabelecendo a ligação com o sul. 


5. DESCRIÇÃO SUMÁRIA DAS MEDIDAS SINTÁCTICAS DOS MAPAS 
ANALISADOS 


As dimensões morfológicas dos mapas axiais resultantes das transforma- 
ções gráficas operadas nas bases topográficas representadas nos mapas 
referidos (Tinoco, Pombalino, Romântico, Ressano Garcia, Estado Novo e 
Actual) estão resumidamente descritas no Quadro 1. 


As variações nas dimensões dos mapas axiais —- número de linhas e nas 
propriedades configuracionais do sistema são devidas fundamentalmente a 
dois factores: à dimensão da população das linhas axiais e às alterações 
ocorridas na morfologia do tecido da cidade. Como foi referido, a dimensão 
do território cartografado varia consideravelmente nos mapas analisados tra- 
duzindo a evolução e expansão urbana. Esta dinâmica é acompanhada por 
alterações morfológicas detectáveis através da estrutura sintáctica e dimen- 
são da configuração. Assim, nas passagens entre os mapas ou plantas |, |l 
e Ill observa-se um decréscimo no número de linhas que compõem os siste- 
mas. Pelo contrário, nos mapas seguintes, dá-se um aumento substancial 
na dimensão do sistema, em particular no que se refere aos dois últimos 
momentos de análise. Os limites da área de estudo sofrem sucessivas 
ampliações até estabilizarem no momento descrito pelo mapa V. Nos mapas 
ou plantas |, Il e Ill, as alterações ocorridas em função da reestruturação do 
tecido e rectificação de traçados traduzem um aumento na acessibilidade 
global. A integração média dos sistemas aumenta de .447 para .754. Nos 
momentos seguintes a acessibilidade global diminui em particular nos dois 
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Fig. 6 Planta Actual 
Planta topográfica de Lisboa, levantada, desenhada e publicada pelo Instituto Geográfico 
do Exército, 1993 


sistemas finais. Este facto relaciona-se com o aumento significativo na 
dimensão do sistema, resultado não só, da ampliação dos limites urbanos 
como do prenchimento de espaços no interior da malha. O padrão de aces- 
sibilidade local dos sistemas, obtido em função do valor do Raio 3, segue o 
padrão acima descrito. As propriedades locais dos sistemas, referenciadas 
aos valores da conectividade e controle, mostram um maior investimento ao 
nivel de relações locais entre os espaços nos sistemas | e Ill em que o coe- 
ficiente de correlação entre as medidas é respectivamente .815 e .801. 


6. ANÁLISE DOS MAPAS AXIAIS PELAS MEDIDAS SINTÁCTICAS 
LOCAIS E GLOBAL 


Com a finalidade de completarmos a análise elaborada anteriormente sobre 
os Mapas Axiais que descrevem, em seis momentos, a evolução morfoló- 
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Planta |  Plantall Planta lll Planta lV Planta V Planta VI 
Tinoco  Pombalino Romântico Ressano Estado Actual 


Garcia Novo 


| n.º de linhas 


INTEGRAÇÃO 


INTELIGIBILIDADE 
GLOBAL 
Int,(Con 
LOCAL-R3,Con 


R3,Cnt 
Con,Cnt 


Quadro | — Valores médios, minimos e máximos da Integração, Raio 3, Conectividade e 
Controle bem como da Inteligibilidade dos mapas de linhas axiais relativos as Plantas 
deTinoco, Pombalino, Romântico, Ressano Garcia, Estado Novo e Actual. 


gica da cidade de Lisboa iremos analisar comparativamente o desempenho 
das medidas sintáticas locais — Conectividade, Controle e Raio 3 — bem 
como global — Integração — destes mapas. 


Na Fig. 13 e no Quadro | podemos constatar a evolução dos valores 
médios, mínimos e máximos destas medidas para as linhas axiais de cada 
mapa. 


O CONECTIVIDADE MIN | [LICONECTIVIDADE MED — ACONECTIVIDADE MAX 


0 1 2 Mod 5 O, 
TINOCO POMBALINO ROMÂNTICO RESSANO ESTADO ACTUAL 
GARCIA NOVO 


O CONTROLE MIN O CONTROLE MED À CONTROLE MAX 


0 1 2 3 4 5 6 7 
TINOCO POMBALINO ROMÂNTICO RESSANO ESTADO ACTUAL 
GARCIA NOVO 


Fig. 13 Evolução dos valores médios, mínimos e máximos da Integração, Raio 3, Conecti- 
vidade e Controle das linhas axiais dos mapas de Tinoco, Pombalino, Romântico, Ressano 
Garcia, Estado Novo e Actual. 


Em relação às medidas locais podemos constatar uma evolução semelhante 
para os mapas analisados. Com efeito, a Conectividade, a medida mais 
local, apresenta valor máximo para o mapa do Estado Novo e menor para o 
mapa de Tinoco o que significa maior número de interseções para aquele e 
menor para este. Os valores médios da Conectividade e Raio 3 são meno- 
res para o mapa de Tinoco, aumentam progressivamente até ao período 
Romântico, diminuindo após a intervenção de Ressano Garcia sobre o 
tecido urbano de Lisboa. 


Morfologicamente esta evolução caracteriza-se pela introdução generalisada 
da grelha como instrumento de desenho urbano na conformação da cidade 
que começa no Bairro Alto, já presente, ainda que parcialmente, no mapa 
de Tinoco. Esta grelha é plenamente assumida no mapa Pombalino e no 
Romântico apresentado-se, posteriormente, como um tipo de intervenção 
mais fragmentada e pontual como ocorre a partir do mapa de Ressano Gar- 
cia. O aumento e declínio destas medidas locais médias pode ser intreper- 
tado também pelo aumento de escala métrica que ocorre a partir do mapa 
do Estado Novo que introduz grandes avenidas quando comparadas com as 
dos mapas anteriores. 


Os valores máximos destas medidas também apresentam consistência evo- 
lutiva pois são menores para os mapas de Tinoco e Pombalino quando 
comparados com os do mapa do Estado Novo mostrando a importância 
local que assume a morfologia urbana neste último mapa. Os grandes eixos 
que conformam a cidade — Passeio Público, Avenidas da Liberdade, Almi- 
rante Reis e República — bem como o aparecimento de bairros morfologica- 
mente homogéneos, como é o caso das Avenidas Novas, contribuem para 
esta situação. 


É, no entanto, em relação à medida global — a Integração — que podemos 
notar maiores variações de comportamento em relação à caracterização 
morfológica destes mapas. Se reparmos na Fig. 13 podemos observar que 
os valores médios e máximos desta medida aumentam, progressivamente, 
desde o mapa de Tinoco até ao mapa de Ressano Garcia apresentado 
depois um acentuado declínio até ao mapa Actual. 


Este resultados terão de ser cuidadosamente interpretados na medida em 
que a Integração é função da distância axial —- ou de mudança de direcção 
da orientação dos percursos mínimos — de todas as linhas axiais em relação 


59 


60 


as restantes. Por outras palavras, a Integração é uma função da dimensão 
do sistema dado pelo número de linhas axiais presentes em cada um deles. 


Com efeito, se repararmos na Fig. 14, onde se apresentam os resultados da 
regressão entre a Integração média de cada mapa axial com a respectiva 


y = -3.466E-5x + .689, R-squared: .357 
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Fig. 14 Regressão da Integração média com o Nº de linhas axiais para os sistemas anali- 
sados com e sem o mapa de Tinoco. 


dimensão, podemos constatar que o coeficiente de determinação é de 0,357 
com 95% de confiância e que o mesmo cresce para 0,771 se lhe retirarmos 
a observação correspondente ao mapa axial que representa o mapa de 
Tinoco. 


Estes resultados mostram que estamos, provavelmente, em face de um 
artefacto na medida em que a Integração é uma função da dimensão de 
cada sistema de linhas axiais analisado. No entanto, estes resultados cha- 
mam-nos também a atenção para a possível ocorrência de dois tipos morfo- 
lógicos distintos de tecido urbano — o medieval com ruas estreitas e tortuo- 
sas com uma conformação quase labiríntica e a que resulta da intervenção 
Pombalina sob a forma de uma grelha reticulada que posteriormente é assu- 
mida como forma preponderante de desenho da cidade mesmo que este se 
manifeste com a configuração de grelha deformada pelas sucessivas “cola- 
gens, cozeduras e acrescentos” que podemos constatar na sua malha. 


Estas considerações clarificam-se (ver Fig. 15) se fizermos uma regressão 
entre a Inteligibilidade de cada mapa axial com a respectiva dimensão. Pode- 


mos, assim, constatar que, neste caso, passamos de um coeficiente de deter- 
minação de 0,303 para 0,948, respectivamente, com e sem o mapa de Tinoco. 
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Fig. 15 Regressão da Inteligibilidade com o Nº de linhas axiais para ossistemas analisados 
com e sem o mapa de Tinoco. 


Por outras palavras, se a Inteligibilidade dos sistemas espaciais diminui à 
medida que nos aproximamos da época actual isso, em grande parte, é 
devido ao aumento das dimensões de cada sistema ao longo do tempo res- 
tando, como hipótese a ser confirmada, se ocorrem ou não dois ou mais 
periodos morfológicos na evolução urbana da cidade de Lisboa que sejam 
dectáveis por aferição das suas propriedades sintáticas para além da sua 
dimensão. 


A evolução do Coeficiente de Variação da Integração (ver Fig. 16) pode for- 
necer-nos uma pista na medida em que o mesmo é invariante em relação à 
população de linhas axiais que compõem cada sistema. 

O Coeficiente de Variação, dado pela divisão entre o desvio padrão e a 
média da Integração compara, de forma padronizada, a dispersão desta 
variável pelos seis mapas independentemente da dimensão de cada popula- 
ção de linhas axiais. 

Assim, podemos observar que este coeficiente diminui de forma significativa 
do mapa de Tinoco para o Pombalino, que aumenta ligeiramente nos mapas 
Romântico e de Ressano Garcia para sofrer novo decréscimo no mapa do 
Estado Novo para voltar a aproximar-se, no mapa Actual, do valor já alcan- 
çado no mapa Romântico. 
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TINOCO  POMBALINO ROMÂNTICO RESSANO ESTADO ACTUAL 
GARCIA NOVO 


Fig. 16 Coeficiente de Variação da Integração para os mapas analisados. 


Podemos confirmar, desta maneira, que a Integração das linhas axiais dos 
mapas de Tinoco e do Estado Novo caracterizam dois períodos morfológi- 
cos da cidade de Lisboa sendo o primeiro descrito como relativamente 
menos uniforme do que o segundo onde a variação da Integração é mínima 
quando comparada com a dos restantes mapas. 


Se a integração média de 0,447 do mapa de Tinoco revela uma cidade algo 
segregada o mapa do Estado Novo apresenta um valor próximo daquele de 
0,490 (ver Quadro |) mas com quase o triplo das linhas axiais que estão 
presentes no primeiro mapa. Isto significa que, em termos comparativos, 
acabamos por ter uma cidade quase tão segregada como a do período 
medieval, muito mais extensa e com menor variabilidade na integração das 
suas linhas axiais. Enquanto no período medieval se promoveu uma maior 
variação morfológica do seu tecido no mapa do Estado Novo este apre- 
senta-se muito menos diferenciadado do que anteriormente. 


Esta constatação mostra-nos que a extensão, melhoramento e consolidação 
do programa de obras públicas encetado durante a vigência do Estado Novo 
até 1948 contribuí para uma homogeneização morfológica do tecido urbano 
só comparável à que ocorreu durante a intervenção Pombalina, após o ter- 
ramoto, mas com uma diferença fundamental: a intervenção do Estado Novo 
requalifica a cidade para limiares de variação morfológica que ficam aquém 
do Plano Pombalino constituindo-se, assim, num dos momentos de cresci- 
mento menos diververgentes, sob o ponto de vista morfológico, de toda a 
história da cidade de Lisboa desde o período medieval. 


7. ANÁLISE DOS MAPAS AXIAIS E DOS SEUS CASCOS DE INTEGRA- 
ÇÃO 


Os mapas axiais de Lisboa, desde o período medieval até ao actual, são 
reveladores das tranformações morfológicas que o tecido urbano da capital 
sofre ao longo de mais de três séculos da sua existência. 


Certamente que o mapa que melhor caracteriza a estrutura urbana da 
cidade ao longo destes períodos, é o relativo à medida de Integração, pois 
relaciona cada linha axial com as restantes fornecendo indicações acerca 
dos espaços urbanos mais integrados e segregados, bem como do respec- 
tivo casco urbano, isto é, do conjunto de linhas axiais que se apresentam 


com valores de integração indênticos ou superiores aos percentis mais ele- 
vados (5%, 10% e 20%). 


O modo como o casco urbano se posiciona em relação à cidade e a forma 
mais ou menos fragmentada ou coesa de como se desenvolve territorial- 
mente são profundamente reveladoras da evolução da morfologia da cidade 
que desde o perídodo medieval, de uma cidade voltada para o rio se trans- 
forma numa cidade que se expande para norte, consolidando-se como sede 
de uma área metropolitana com um centro direcional formado pelas Aveni- 
das Novas situadas no prolongamento e expansão das malhas urbanas mais 
antigas que se desenvolvem à volta das Avenidas Almirante Reis e da Liber- 
dade. 


Se observarmos o Mapa de José Luis Tinoco bem como o correspondente 
mapa axial em relação à medida de integração (ver Figs. 1 pág. 55 e 7) po- 


Fig. 7 Mapa axial da planta de Tinoco relativo à medida de Integração 


demos constatar que o conjunto de linhas axiais mais integradas, represen- 
tadas neste último mapa a traço mais escuro, dispõem-se e concentram-se, 
de forma coesa, paralelamente ao rio fazendo uma única penetração para o 
interior pela Rua Nova dos Ferros. Isto mostra a importância morfológica 
atribuída às relações da cidade com o rio bem como ao comércio represen- 
tado por aquela Rua Nova apesar de, como nos lembra Oliveira Marques 
(1994) a cidade já ter anteriomente ao Séc. XVII “quatro ou cinco centros ou 
pólos aglutinadores: a alcáçova, ainda centro militar e político, onde vivia o 
rei e se situava a corte; a catedral ou Sé, mantendo as suas funções religio- 
sas ímpares, embora ameaçada pela proliferação de mosteiros e conventos 
mais periféricos; a Ribeira, o Rossio; e a Rua Nova na Baixa”. Deste con- 
junto de pólos aglutinadores apenas a Ribeira e a Rua Nova pertencem ao 
casco urbano das linhas axiais mais integradas ficando os restantes em 
linhas mais segregadas, representadas no respectivo mapa axial pelas tona- 
lidades mais claras confirmando que Lisboa dos meados Séc. XVII é uma 
cidade morfologicamente organizadaem função do Rio. O mapa evidencia 
ainda a importância do Paço Real situado no Terreiro do Paço que, conjun- 
tamente com a Rua Nova, as Tercenas e os estaleiros de navios, contri- 
buem para que o centro da cidade se defina na zona mais ribeirinha da 
Baixa, confirmado pelo casco urbano do mapa axial que se apresenta coeso 
e integrado. O Rossio, situado mais a norte, ocupa uma posição periférica 
em relação ao casco urbano ao contrário do Terreiro do Paço que se loca- 
liza na confluência das linhas de maior integração. A morfologia da organiza- 
ção do espaço urbano neste período consegue reproduzir as relações soci- 
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ais que se operam na sociedade, na medida em que o espaço urbano se 
revela como uma instância de produção e reprodução dessas mesmas rela- 
ções. Com efeito, os espaços urbanos mais integrados correspondem a um 
maior incremento e intensificação das actividades urbanas que ocorrem ao 
longo deste período. A Lisboa oriunda da época medieval representada no 
mapa de Tinoco, caracteriza um período morfológico da evolução urbana de 
Lisboa, profundamente distinto dos que se lhe sucederam. 


É com as consequências do terramoto de 1755 que se opera uma ruptura 
morfológica na organização espacial da cidade pela adopção do Plano Pom- 
balino de profundas implicações até aos dias de hoje. 


Não se trata somente de produzir um novo plano que rectifique e alargue os 
arruamentos estreitos e tortuosos que foram legados pelo período anterior, 
mas de propor um novo desenho para a cidade e, em particular, para o seu 
centro urbano. 


Se as questões de ordem que o novo plano suscita podem ser encaradas 
como uma revitalização da Lisboa que sai dos escombros pelo respeito à 
localização das permanências, principalmente do Rossio e do Terreiro do 
Paço que se mantêm praticamente inalterados, o mesmo já não se passa 
com a estrutura urbana que emerge do novo Plano. 


Ao analisarmos a Planta Pombalina bem como o correspondente mapa axial 
relativo à medida de integração (ver Figs. 2 pág. 56 e 8) podemos constatar 


Fig. 8 Mapa axial da planta Pombalina relativo à medida de Integração 


que o casco urbano é distinto do anterior. Não só se expande para o inte- 
rior, deixando adivinhar a importância que o Passeio Público passará a 
assumir, como promove a ligação entre a Lisboa Oriental — da zona segre- 
gada do Castelo — com a Ocidental — da zona do Bairro Alto. Trata-se de 
uma nova malha reticular que promove a integração das diferentes partes 
da cidade. Do Bairro Alto pela saída de Santos e pela Rua do Alecrim e da 
zona do Castelo pela Rua da Madalena, afastando-se do rio e quase que 
“voltando as costas” a este. 


As relações intra-urbanas passam a apresentar um significado não assu- 
mido pela Lisboa oriunda do período Medieval. O Terreiro do Paço continua 
a existir como espaço simbólico e de representação do poder do estado 
mas a sua importância na estrutura urbana é menor, já não faz parte do 
casco mais integrado e situa-se periféricamente em relação a este apesar 
de ser intersectado por linhas axiais oriundas do Rossio que cobrem as 
Ruas do Ouro, Augusta e da Prata. É uma nova cidade que promove, pela 
sua morfologia espacial, as relações intra-urbanas e, de algum modo, se 


fecha em relação ao rio apesar de o Terreiro de Paço localmente promover 
essa abertura. O dominio da frente de rio sobre a cidade opera-se local- 
mente mas não participa do casco urbano mais integrado. 

O Rossio deixa de ter uma posição periférica para passar a fazer parte do 
casco urbano da cidade constituindo-se como pólo alternativo ao Terreiro do 
Paço. À separação das escalas monumental e gregária, de representação 
do estado e de actividades quotidianas, são assumidas agora, de forma 
separada, pelas praças do Terreiro do Paço e pelo Rossio. 


O Passeio Público, desenhado sobre terrenos de quintas, não é incorporado 
no casco urbano mais integrado mas situa-se adjacente a este sugerindo o 
prolongamento da Baixa Pombalina ao longo de um vale — o Valverde. Se 
bem que a escala gregária de encontro e lazer do empreendimento definisse 
o destino do novo Passeio, a sua organização como espaço alinhado de 
árvores com a “aparência de quinta no fundo da província ou larga alameda, 
assaz monacal” (França 1990) não favorecia a sua utilização pelos lisboetas. 
A sub-utilização deste espaço público não é estranha a sua posição perifé- 
rica em relação ao casco urbano mais integrado do Plano Pombalino. 

É na Lisboa Romântica (ver Figs. 3 pág. 56 e 9) que podemos observar que 
o casco urbano já comparece constituído pelo Passeio Público dando “cate- 


Fig. 9 Mapa axial da Planta Romântica relativo à medida de Integração 


goria social e urbana ao sítio, justificando edificações e valorizando palácios 
pelos alinhamentos” (J. R. Ferreira, 1994) participando activamente na orga- 
nização morfológica da cidade que se expande para norte. Com efeito, o 
próprio Passeio organiza-se morfológicamente num anel de alinhas axiais 
mais integradas que o envolvem contrabalançando a existência de outra 
estrutura axial anelada formada por linhas axiais pertencentes à Baixa Pom- 
balina, Bairro Alto, saída de Santos e Rua do Alecrim. 


Estas estruturas aneladas de linhas axiais integradas reproduzem aquilo que 
Hillier e Hanson (1984) designam de “roda deformada”, isto é, do conjunto 
de linhas axiais que se organizam como radiais de uma configuração, apro- 
ximadamente, em forma de roda. Segundo estes autores a roda deformada 
é característica dos assentamentos urbanos mais expontâneos, ou seja, 
daqueles que estão sujeitos, predominantemente, a regras de organização 
local e não global — à escala de todo o assentamento — como é o caso da 
Baixa Pombalina. 


Nos mapas da Lisboa Romântica (ver Figs. 3 e 9) podemos constatar a for- 


mação dessas estruturas em forma de roda deformada indicando as mes- 
mas a direcção da futura expansão urbana, mais orgânica e local, para 
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norte, e ao longo do rio, para ocidente. É, com efeito, neste período que se 
implantam o jardim da Estrela, o miradouro de S. Pedro de Alcântara e do 
Príncipe Real, confirmando, os dois ultimos, um percurso urbano do Cais 
Sodré ao Rato. Igualmente surgem a urbanização do Aterro da Boavista ao 
longo da faixa ribeirinha que será a futura 24 de Julho. 


Estas formas aneladas do casco sugerem um desenvolvimento urbano mais 
diversificado quando comparado com o período precedente. Enquanto a Baixa 
surge a partir dos planos dos engenheiros militares Pombalinos, desenhada 
pelo poder centralizador do estado — onde predominam as regras globais 
sobre as locais para a construção da nova cidade — as novas extensões do 
casco urbano ao longo do Passeio Público e do anel de linhas axiais Integra- 
das que passam pela Praça do Cais do Sodré indicam um tipo de crescimento 
sujeito a regras locais geradoras da expansão da cidade Romântica — mais 
permeáveis a investimentos pontuais e privados, promovendo uma maior 
diversidade de intervenções arquitectónicas. 


Estas tendências acentuam-se de alguma maneira no Plano de Ressano Gar- 
cia (ver Figs. 4 pág. 57 e 10) na medida em que a estrutura anelada de linhas 
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Fig. 10 Mapa axial da Planta Ressano Garcia relativo à medida de Integração 


axiais mais integradas que passam pelo Cais do Sodré se mantém relativa- 
mente intacta. A grande diferença ocorre pela exclusão do que outrora foi o 
Passeio Publico e a sua substituição pela Avenida da Liberdade configurando- 
se a linha axial que a cobre como pertencente ao casco urbano da cidade. 
O anel que envolvia aquele Passeio desaparece para dar lugar a uma longa 
linha axial que expande o casco urbano para norte. isto sugere tratar-se de 
um compromisso com os dois momentos anteriores : por uma lado, a preser- 
vação de uma roda deformada que se desenvolve ao longo do rio para oci- 
dente e, por outro, o apontar de uma nova direccionalidade que não estava 
presente no desenho de Lisboa. O germe do futuro centro direcional, constitu- 
ido por um casco urbano de linhas axiais bem integradas, já está presente no 
plano de Ressano Garcia para Lisboa. Esta confirmação contraria a ideia colo- 
cada por V.M. Ferreira (1994) de que a nova Avenida seria uma referência 
urbanistica ao “boulevard Haussmaniano” e não um modelo de revovação 
urbana como propõem M.J. Rodrigues (1979) e J. M. Fernandes (1980). Trata- 
se, com efeito, da génese da nova cidade que se irá expandir a partir das 
Avenidas Novas embora ainda constituam áreas segregadas no Plano de Res- 
sano Garcia. Se bem que o desenho da nova malha destas Avenidas retome 
o reticulado geométrico da Baixa Pombalina e, portanto, traduza um modelo 


gerado por leis globais de desenho e ordenamento urbano as funções 
desempenhadas por ambas são, ainda, distintas. As primeiras como espaço 
de expansão urbana espacialmente segregado da escala monumental da 
cidade e a segunda como intervenção centralizadora e conformadora das 
relações espaciais da sociedade civil com o novo estado absolutista. 


É na simbiose destes “momentos” de desenho urbano que se consolida a 
participação do Estado Novo (ver Figs. 5 pág. 58 e 11) pelo retomar do Plano 
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Fig. 11 Mapa axial da planta do Estado Novo relativo à medida de Integração 


das Avenidas Novas pelo Eng. Duarte Pacheco nos finais dos anos 20 com 
a implantação do complexo do IST. O fomento das Obras Públicas associ- 
ado a uma política administrativa autoritária que permitiu ao Município de 
Lisboa apropriar-se, entre 1935 a 1949, de cerca de um terço da área total 
concelhia conduziu à formação de uma enorme bolsa de terrenos que serviu 
de base para o lançamento de novos Bairros e equipamentos bem como 
pela consolidação de situações expectantes. 


É no Estado Novo que o casco urbano de Lisboa se “desloca” para norte 
englobando as linhas axiais que cobrem, por um lado, as Avenidas da Liber- 
dade, António Augusto Aguiar, Fontes Pereira de Melo e República e, por 
outro, a Almirante Reis, a Avenida de Roma e a Gago Coutinho sendo 
ambas unidas por transversais que passam pela Avenida João XXI, Avenida 
Estados Unidos da América e Morais Soares. 


Todo o casco urbano de Lisboa se reconstitui e expande para norte dei- 
xando a Baixa de apresentar a centralidade que anteriormente tivera bem 
como a área envolvente ao Cais do Sodré deixará de se configurar também 
como pertencente áquele casco. É uma nova centralidade que leva a cidade 
para norte, propondo duas grandes linhas de penetração que se conformam 
aos vales, perpendiculares ao rio não deixando, contudo, as linhas axiais 
que constituem este casco de formarem aneis. 


Não se trata já de uma configuração do casco urbano em roda deformada, 
mas de uma organização em múltiplos aneis de linhas axiais integradas que 
proporcionam percursos e usos alternativos formando o centro direcional de 
uma nova metrópole no qual se implantarão serviços e comércio com uma 
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lógica social de ocupação do espaço que confirma o novo casco urbano da 
cidade a norte (veja-se Gaspar, 1976). Além disso, a escala metrica deste 
casco altera profundamente as relações anteriormente estabelecidas pelo 
aumento da dimensão das linhas axiais. Compare-se a dimensão linear das 
Ruas do Ouro, Augusta e da Prata com cerca de 550 metros de compri- 
mento com as da Almirante Reis com 3 100 metros e a Avenida da Repu- 
blica com 1 700 metros. 


Outra escala impera neste novo casco muito mais voltado para a definição 
de uma direcionalidade norte do que para a concentração das actividades 
urbanas mais intensas numa área restrita. Trata-se de uma nova cidade 
assumida morfológicamente como centro direcional constituido por aneis 
múltiplos de linhas axiais integradas com uma escala métrica nunca antes 
desenhada onde predomina uma morfologia que “agarra e coze” as diferen- 
tes partes do tecido urbano mais periférico e segregado. 


No periodo actual (ver Figs. 6 pág. 59 e 12) podemos constatar o reforço 
das tendências anteriormente detectadas no Plano do Estado Novo através 


Fig. 12 Mapa axial do momento actual relativo à medida de Integração 


da clarificação do casco urbano formado pelas linhas axiais que cobrem as 
Avenidas Almirante Reis e Gago Coutinho, por um lado, e a da Liberdade, 
Fontes Pereira de Melo e da República, por outro. As linhas axiais transver- 
sais a estas últimas deixam de configurar anéis múltiplos simplificando a 
conformação do casco urbano que se apresenta agora com uma clara direc- 
cionalidade norte. Este sentido é, em parte consequência do preenchimento, 
por novos empreendimentos, dos “hiatos” urbanos deixados no Mapa do 
Estado Novo: Olivais Norte e Sul, Chelas e Telheiras propiciam um novo 
casco urbano definido de uma forma ramificada para norte configurando-se 
a Av. João XXI como a única via urbana coberta por linhas axiais integradas 
que se dispõem transversalmente ao desenvolvimento norte do casco 
urbano. A Baixa Pombalina deixa de ser definitivamente o centro morfoló- 
gico de Lisboa para passar a ser mais um dos seus bairros. 


Em resumo, desde o Plano de Tinoco, em que a cidade se voltava intencio- 
nalmente para o rio, assistimos a uma progressiva orientação do casco 


urbano para norte configurando-se este, nos dias actuais, como centro dire- 
cional plenamente assumido. 


Qualquer novo plano que requalifique as relações da cidade com rio terá de 
levar em consideração a nova morfologia da cidade que se foi estruturando 
após o terramoto de 1755 e entender que somente poderá aspirar a uma 
nova conformação se propuser alterações locais que não interfiram global- 
mente na morfologia urbana já consolidada da cidade de Lisboa. 


CONSIDERAÇÕES FINAIS 


Como resultado das análises elaboradas podemos detectar a ocorrência de 
dois períodos morfológicos na cidade de Lisboa: antes e após o terramoto. 
Podemos mesmo dizer que estamos em presença de duas cidades que se 
caracterizam de modo diversificado em relação à dimensão global das suas 
linhas axiais bem como à configuração e posição do casco urbano em rela- 
ção a todo o seu tecido urbano. 


A mais antiga, mais segregada mas menos homogénea em termos de Inte- 
gração das suas linhas axiais, e a mais moderna, menos segregada após o 
terramoto mas progressivamente menos integrada, culmina numa menor 
variabilidade do seu tecido urbano ao qualificar o seu território de forma 
relativamente indiferenciada, principalmente no período do Estado Novo. 


Conclui-se pelo momento notável da alteração da morfologia da cidade dada 
pela intervenção Pombalina, coincidente com a modificação operada na 
sociedade portuguesa, só possivel devido à existência de uma catástrofe 
natural aliada a um forte domínio do poder centralizador do Estado. Esta 
alteração é contudo acompanhada de uma diminuição significativa da inteli- 
gibilidade espacial do tecido urbano pelas intervenções urbanas que suces- 
sivamente ocorreram desde a catástrofe de 1755, principalmente pelas 
modificações introduzidas no período do Estado Novo que promovem uma 
acentuada homogeneização morfológica do seu tecido urbano bem como a 
consolidação de um centro direcional que se afasta progressivamento do rio. 


REFERÊNCIAS BIBLIOGRÁFICAS 


Fernandes, J. M. (1980) “Alguns apontamentos sobre urbanismo alfacinha. Da Baixa 
ao Campo Grande, do Séc. XVIII ao Séc. XX". Arquitectura, ano 11, 4.º Série, 
138. 


Ferreira, J.R. (1994) “Passeio Público”. Verbete do Dicionário da História de Lisboa, 
Direcção de F. Santana e E. Sucena. Carlos Quintas & Associados — Consultores 
Ltda..Lisboa. 


França, J. A. (1976) “Lisboa Oitocentista”. Academia Nacional de Belas Artes/ Funda- 
ção Calouste Gulbenkian. Exposição Documental. Lisboa. 


França, J. A. (1987) “Lisboa Pombalina e o Iluminismo”. Bertrand Editora. Lisboa. 


França, J. A.(1990) “ A Arte em Portugal no Século XIX”, Lisboa, vol.ll, 2.º Ed., Ber- 
trandEditora 


Gaspar, J. (1976) “A Dinâmica do Centro Funcional de Lisboa”. Livros Horizonte. Lis- 
boa. 


Hillier, B. e Hanson, J. (1984) “The Social Logic of Space”. Cambridge University 
Press. Cambridge.U.K.. 


Krúger, M.J.T. (1990) “On Node and Axial Grid Maps: Distance Measures and Rela- 
ted Topics"Working Paper. The Bartlett School of Architecture. University College 
Lodon. 


Marques, A.H.O. (1984) “Lisboa Evolução; Séculos XIl a XV”. Verbete do Dicionário 
da História de Lisboa. Direcção de F. Santana e E. Sucena. Carlos Quintas & 
Associados — Consultores Ltda..Lisboa. 


69 


70 


Moita, |. (1982) “Lisboa e o Marquês de Pombal”, vol.ll, Lisboa, Museu da Cidade, 
CML 

Moita, |. (1983) “Lisboa Quinhentista — A Imagem e a Vida da Cidade”. Direcção dos 
Serviços Culturais da Câmara Municipal de Lisboa. Museu da Cidade. Lisboa. 

Nunes da Silva, C. (1987) “Planeamento Municipal e Organização do Espaço em Lis- 
boa: 1926-1674', Tese de Mestrado CEG-UL, Lisboa 

Rodrigues, M.J. (1979) “Tradição, transição e mudança. A produção do espaço 
urbano na Lisboa oitocentista". Boletim da Assembleia Distrital. Lisboa. 

Santana, F. (1994) “Lisboa Evolução: Periodo da Expansão”. Verbete do Dicionário 
da História de Lisboa, Direcção de F. Santana e E. Sucena. Carlos Quintas & 
Associados — Consultores Ltda..Lisboa. 


ÍNDICES ANTERIORES 


N.º 1/94 - Junho 1994 


ARTIGO CONVIDADO 


Aceleradores de Partículas por Interacção Laser - Plasma: J.T. Men- 
donça e L. Oliveira e Silva 


COMUNICAÇÕES CIENTÍFICAS 


Termodinâmica para Engenharia Mecânica: Finalidades e Objectivos: 
João Ventura 


Chemistry, Geology and Materials Science of Geothermal Systems: 
M.O.S.P. Caldeira and C.A.C. Sequeira 


Nil Contribution of the Cosmological Constant to the Plaine Gravita- 
tional Waves in General Relativity: K.S. Adhav and T.M. Karade 


A Cilindrically Symmetric Cosmological Model in Bimetric Relativity: 
Shailendra Deo 


A Reflectmetric System for the Tokamak Isttok: P. Varela, A. Silva, S. 
Verganota, F. Serra, M.E. Manso, J.A. Cabral and C.A.F. Varandas 


71 


The Difficulty of Attributing a Physical Significance to Quantities 
Work and Heat in the First Principle of Thermodynamics. The 
Resolution of an Entropy Maximization Controversy: Rodrigo de 
Abreu 


Detecção de Solutos em Cromatografia lónica - Detector Condutimé- 
trico: J. Costa Pessoa 


Natureza e Estrutura da Revista Técnica 


N.º 2/94 - Setembro 1994 


INTERVENÇÕES NA CERIMÓNIA COMEMORATIVA A DUARTE PACHECO: 
Professor Diamantino Freitas Gomes Durão 
Presidente do IST 


Professor Jorge Borges de Macedo 
Director do Arquivo Nacional/Torre do Tombo 


Professor E.R. de Arantes e Oliveira 
Director do LNEC 


Professor A. Simões Lopes 


Reitor da Universidade Técnica 


Engenheiro António Silva Teixeira 
Presidente da Direcção da AAAIST 


Engenheiro Ferreira do Amaral 


Ministro das Obras Públicas, Transportes e Comunicações 


TESTEMUNHOS SOBRE VIDA E OBRA DE DUARTE PACHECO: 


Duarte Pacheco um Homem do futuro: Luís Guimarães Lobato 


Posicionar o Instituto Superior Técnico para o início do Século XXI: 
P. Conceição, D.F.G. Durão e M.V. Heitor 


Duarte Pacheco - Acção, Engenharia, Elevação: M. de Abreu Faro 


Concepção Arquitectónica de Instituto Superior Técnico: Madalena 
Cunha Matos (Arquitecta - Departamento de Engenharia Civil do IST, 
Universidade Técnica de Lisboa) 


N.º 3/94 - Dezembro 1994 


ARTIGO CONVIDADO 


A Fusão Nuclear Controlada, o Programa Comunitário da Fusão e a 
Participação Portuguesa: J.A.C. Cabral, M.E. Manso e C.A.F. Varandas 
72 | | 
COMUNICAÇÕES CIENTÍFICAS 


The Ideal Gas Equation P = ou and the Energy-Entropy Principle: 
Rodrigo de Abreu 


Electrochemical Properties of Polymer Electrolytes: L.P.S. Araújo and 
C.A.C. Sequeira 


Redes Neuronais em Robótica: Pedro Quintas Aguiar e Luís Borges de 
Almeida 


The Concepts of “Gorce” and “Gheat”: What is the Problem with the 
First Law of Thermodynamics?: Rodrigo de Abreu 


Os Sistemas de Vácuo de Limpeza da Câmara e de Injecção de gás 
do Tokamak ISTTOK: C.A.F. Varandas, C. Freitas, V. Prego, A. Soa- 
res, B. de Groot, F.M. Dias, B.B. Carvalho, J. Sousa, H. Fernandes, 
S. Magalhaes e J.A.C. Cabral 


Térmica de Edifícios: uma Solução Integrada para Distintos Utiliza- 
dores: J.M.S. Dionísio 


Detecção de Solutos em Cromatografia lónica - Detector Espectrofo- 
tométrico: J. Costa Pessoa 


N.º 1/95 - Março 1995 


ARTIGO CONVIDADO 


A Ciência dos Sistemas Complexos: Rui Dilão 


CIÊNCIA E TECNOLOGIA: NOTAS BREVES 


Heinrich Hertz: M. de Abreu Faro 
Prémio da Boa Esperança 1993 atribuido a professor do IST 
Plasmas Coloidais: David P. Resendes 


Novos Materiais para aplicações em Engenharia: Rogério Colaço 


COMUNICAÇÕES CIENTÍFICAS 


Salinidade Atmosférica: um vector natural de risco: A.M.G. Pacheco, 
A.J. Sousa, R. Figueira, C. Branquinho, F. Catarino 


Thin Film Solid state microbatteries: L.P.S. Araújo, M.J.C. Plancha and 
C.A.C. Sequeira 


Detecção de Solutos em Cromatografia lónica - Detecção Espectro- 
fotométrica Indirecta: J. Costa Pessoa 


A Mecânica Quântica de Heisenberg: João Paulo Neves Monteiro dos 
Santos 


Equação de Duffing e Método de Liapunov-Schmidt: Diogo Luís 
Aguiar Gomes 


N.º 2 / 95 - Junho 1995 


ARTIGO CONVIDADO 


Funções Delta no Universo da Análise Não-Standard: J. Sousa Pinto 


NOTAS BREVES 


Estudo geofísico e vigilância de vulcões activos - o vulcão do fogo: 
Sandra Heleno da Silva 


Teoria de Perturbações e ressonância - o pêndulo simples: João 
Paulo Santos 


A utilização de polímeros em fotónica: Luís Miguel Lapão 


COMUNICAÇÕES CIENTÍFICAS 


Salinidade Atmosférica: Primeiras etapas de um programa de biomo- 
nitorização: R. Figueira, C. Branquinho, A.M.G. Pacheco, A.J. Sousa, 
F. Catarino 


73 


A Produção Industrial de Anidrido Ftálico a partir de O-Xileno sobre 
Catalizadores de V,O,/TiO,: C.R. Dias, M. Farinha Portela, J. 
Lemos Pinto 


Sobre a estabilidade de populações iónicas não girotrópicas: A.L. 
Brinca 


Simulação da operação de compressores por programas de cálculo 
automático: João F. P. Gomes 


On the theory of the simple rotating disc electrode: P.S.D. Brito, 
R.M.M. Antunes, C.A.C. Sequeira 


N.º 3/95 - Outubro 1995 


ARTIGO CONVIDADO 


Que Futuro para a Investigação em Engenharia Civil: E.R. Arantes e 
Oliveira 


NOTAS BREVES 


Desenvolvimento de um Inventário das Águas Subterrâneas de Por- 
tugal: Doutor J.P. Cárcomo Lobo Ferreira 


RPM - Rede Inter-Universitária Portuguesa para a promoção e difu- 
são das tecnologias Microelectrónicas: Prof. Epifâneo da Franca e 
Eng.º Rosário Pinto 


74 — COMUNICAÇÕES CIENTÍFICAS 


The Relevance of Non-imaging optics for solar energy: a review: 
Manuel Collares Pereira 


Geometrical relations for non-inertial observers: M. Hermínia Marçal 


Propagação de erros em problemas de mínimos quadrados: Diogo 
Luis Aguiar Gomes 


On the theory of the rotating ring-disc Electrode: P.S.D. Brito, R.M.M. 
Antunes and C.A.C. Sequeira 


Incêndios em núcleos urbanos antigos. Verificação da segurança 
contra incêndios na Mouraria: José Manuel Felizardo Gonçalves 


N.º 4 / 95 - Dezembro 1995 
(Edição comemorativa dos 70 anos da Revista Técnica) 


ARTIGO CONVIDADO 


Ad Multos Annos: J. Fernando de Souza 


COMUNICAÇÕES CIENTÍFICAS 


Modernas concepções da Mecânica: A. Mira Fernandes 


Estructura da Luz - Hipótese dos Quanta, Métodos de Estatística: 
António da Silveira 


Um Vulcão perto de Lisboa entre o Morganhal e Laveiras: Prof. 
Ernest Fleury 


Investigação Científica e Indústria: A. Herculano de Carvalho 


Notas sobre a determinação do Manganês nos seus minérios: 
A. Herculano de Carvalho 


Novo método de cálculo e acerto dos diagramas no projecto das tur- 
binas mixtas (Francis): Manuel P. R. da Costa de Barros 


Estudo sobre o electrodinamómetro de indução: C. Ferrer Moncada 


Movimento permanente gradualmente variado, estudo geral do 
Regolfo: Alberto Abecasis Manzanares 


A Fotoelasticidade no ensino e na técnica: Manuel Mendes da Rocha 


15 


TÉCNICA 


| Revista do Engenharia 
| Public 
de im nara fala o desenvolvimento 


UMA REVISTA DE ENGENHARIA 


A Técnica, Revista de Engenharia, tal como a Associação dos Estudos do Institu. 
Superior Técnico, onde se integra, acompanha, desde os primeiros tempos, 
vida do Instituto. 

Das suas páginas constam artigos científicos e técnicos de flagrante actualidad 
quantas vezes pioneiros, entre nós, de ciência nova e de novas tecnologias. | 

A revista tem uma periodicidade de quatro números por ano e está estruturada e 
perfeita adequação à intensa actividade de investigação e desenvolvimento que hoje s 
processa no campus de Instituto Superior Técnico. 

Compreende, assim artigos convidados, comunicações científicas originais e notícic 
sucintas e rigorosas sobre ciência, tecnologia e desenvolvimento. 

Trata-se, pois, de uma revista cientifica de nível o que não impede a colaboração c 
estudantes desde que devidamente credenciados por investigadores qualificados e c 
reconhecido mérito. 


Torne-se assinante da Técnica. 


Para mais informações contacte-nos pelo telefone 351-1-848 1018 / 8417362 ou fax 351-1-847 40 14 / 849 924 
ou escreva-nos para: Técnica, Revista de Engenharia, Av. Rovisco Pais n.º 1, P-1098-Lisboa CODEX | 
Publicidade: DOSSIER —- Comunicação e Imagem, Lda., R. Eça de Queirós, 18 - 2.º 1000 Lisboa 
Telefs. (01) 355 92 94/5 - 3537321 — Fax/Telef. (01) 3537312 


Reprodução da Maiadouro 


Obra de Tiziano 


Iguais em Tudo, 
Menos na Moldura. 


A qualidade de impressão da Gráfica Maiadouro impressiona mesmo os mais exigentes. Muito 
dificilmente Tiziano distinguiria a sua obra-prima da nossa obra gémea. À prova que a mais avançada 
tecnologia e a completa especialização dos nossos impressores só têm rival na peça original. 
São mais de trinta anos a imprimir qualidade a tudo o que é Livros e Revistas, Catálogos, Desdobráveis, 
Brochuras ou Cartazes, sempre na linha da frente das Artes Gráficas em Portugal. Controlando todos os | 
processos de produção, da fotocomposição e fotografia à montagem, reprodução e acabamento, com o 
rigor e o conhecimento dos grandes mestres. 
É por isso que a assinatura Maiadouro é a melhor garantia de um trabalho perfeito. 
Mesmo sem moldura. 


LAR AFICA 
AAIADOURO, S.A. 


IMPRIME-LHE QUALIDADE. 


SEDE: DELEGAÇÃO: 
Rua Padre Luis Campos, éBé - VERMOIM Avenida Elias Garcia, 162 - 2.º Esg. 
Apartado 1006 - 4471 MAIA Codex 1050 LISBOA 
Tel,: (02] 944 B2 45/84 45 - 948 BB 30/60 04 Tel.: (01) 795 53 00 - 796 92 10 


Fax: [02) 9248 77 07 Fox: [01] 796 92 79 


Ciência Programa Praxis XXI* 


Programa Stride* 
: Programa Interreg I* 
lecnologia apa: 
Investigação Científica e 
p ara Tecnológica (PBICT) 
antecipar Programas Específicos: 
Ciências da Saúde 


o futuro. Preá 


Ciências Sociais e 


Humanas 
Lusitânia 


Ordenamento do 

Território e | 
Desenvolvimento 
Urbano | 


Cooperação 
Universidade- 


Empresa 


Formação e 


Mobilidade de | 


Recursos 


Humanos 


Fundo de Apoio 
à Comunidade 
Científica 
(FACC) 


Cooperação 
Internacional 


em C&T 


Acções de Apoio 


Geral ao Sistema 


g 
E 
É 


Científico 


junta nacional de investigação científica e tecnológica 


Av. D. Carlos |, 126 

1200 LISBOA 

Telets.: 397 9021/7 

Telefax: 60 74 81 

a End. Teleg.: JUNIC-LISBOA 
MINISTERIO DA CIENCIA E DA TECNOLOGIA E-mail : jnictGbgeral, jnict.mailpac.pt 


